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ABSTRACT 

Linear optimal control techniques are used as a synthesis tool to con- 
ceptually design control systems for a large, highly flexible launch vehicle. 
Quadratic performance indices are specified, some of which include a drift 
minimum model, that yield a realizable optimal control law providing ade- 
quate damping of two structural modes as well as satisfactory closed-loop 
speed of response. These designs are shown to be relatively insensitive to 
bending mode shapes and slopes. A theory of optimal control is developed 
for systems possessing parameters whose values are known only on a sample 
space and an optimal compensation network is designed for a launch vehicle 
with an uncertain first bending mode slope. In addition, a theory is devel- 
oped for a design procedure that limits the feedback gains from one or more 
state variables. An experiment is described that obtains optimal solutions 
for systems having one or more feedback gains arbitrarily set to preselected 
value s . 
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SECTION 1 

INTRODUCTION AND SUMMARY 


The design of control systems for launch vehicles must account for 
the effects of structural flexibility and slosh. These effects present more 
and more significant complications as launch vehicles increase in size and 
slenderness. The approach of removing flexibility effects by filtering is 
no longer sufficient, and active control of structural modes is becoming a 
necessity. Thus, a broader demand is placed on the control system de- 
signer. 


The use of linear optimal control theory to meet this demand is the 
subject of this report. The theory has been developed to the point where 
it offers a systematic approach to complex problems of automatic control 
system design. The aim of the study was to show that this theory can be 
used as a practical and effective tool for the design of control systems for 
flexible launch vehicles. 

Linear optimal control theory provides a synthesis procedure for 
linear systems which determines in any particular case a unique closed- 
loop characteristic polynomial. A unique linear combination of the state 
variables is fed back to the control inputs of the plant in accordance with t 
this determination. The theory recognizes the desirability of keeping errors 
in system outputs small, while at the same time using amplitudes of control 
motions that are no larger than necessary. This dual objective is expressed 
in terms of minimizing a performance index, which includes weighted meas- 
ures of the output motions and the control input motions. The approach can 
be explained in terms of the problem of returning a system to equilibrium 
from a disturbance state that is represented by a set of initial conditions. 

The following quadratic performance index is used 

ao 

2 V = J (y'Qy + u.'12u.)dt 
0 

In this expression, y is a vector representing the deviation of the system 
state from the equilibrium state; u is a vector representing the control 
input quantities; and Q and IS are weighting matrices determining the rela- 
tive importance attached to minimizing the variations of y and u respec- 
tively. The variables y and Li are related as follows. The linear constant 
coefficient equations of motion of the vehicle are put in the form 

% = F%+ (3 a 

where % is a vector defining the state of the system, and f and Q are ma- 
trices. The output quantities of interest are determined by the equation 

y 

The problem is to find a 0 > the time history of it which minimizes the per- 
formance index. Then the theory shows that it is possible to generate 
from linear combinations of the state variables oc . This means that it is 
possible to find a feedback matrix ^ such that optimal performance will be 
obtained if we use the control law 

u 0 - -16%. 
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A family of optimal solutions can be found, depending on the choices 
of the weighting matrices Q and 13 . A major advantage of this approach is 
that the theory guarantees that all members of this family will be stable. It 
is found also that the optimal system configurations generally have smooth, 
well-behaved responses. The matter of choosing among the optimal solu- 
tions is a subject of discussion in the body of this report* 

The unique problems associated with the control of a large, flexible 
launch vehicle make the optimal control approach attractive. To date, most 
of the launch vehicle control system design techniques have used filtering of 
one kind or another, such as notch filters, to attenuate the bending mode sig- 
nals in the feedback control loop. This procedure has the effect of gain- or 
phase- stabilizing the poles associated with the bending modes. Although the 
poles are not destabilized, the damping ratios of the bending mode are fre- 
quently not much improved over their open-loop values. Wind or gust dis- 
turbances can still excite the lightly damped modes of the structure. 

As launch vehicles become larger and more slender, the principal 
bending mode frequencies become lower, approaching the desired closed- 
loop natural frequency of the poles originally associated with the rigid-body 
motions. Under these circumstances, it is no longer practical to filter or 
attenuate the bending modes. Instead, it becomes necessary to exert posi- 
tive control over the bending motions of the vehicle. 

Optimal control techniques yield a positive, direct design procedure 
for controlling the bending motions of a large, flexible launch vehicle, while 
at the same time accomplishing the task of controlling the overall rigid-body 
motions of the vehicle. As noted above, a linear optimal control system is 
guaranteed stable despite the complexities introduced by the structural 
modes. More than the assurance of stability is gained by the use of linear 
optimal control techniques. Inherent in the quadratic nature of the perfor- 
mance index is the tendency to penalize large motions of the error and the 
control, yet place little penalty on small errors or deviations. In his book, 

S. S. L. Chang (Reference 4) demonstrates that for the same speed of re- 
sponse, an optimal design requires lower amplitudes of control motion than 
does a conventionally designed system. 

In addition to these advantages, the demands of the quadratic perfor- 
mance index require that closed-loop poles of the system approximate a 
Butterworth filter distribution. The significance of this is that the optimal 
designs provide substantial increases in the damping ratios of the bending 
modes. The elastic motions of the vehicle are blended with the rigid-body 
motions in such a way that the net response of the vehicle tends to be smooth 
and well-behaved to either a command input or other excitations. 

The application of linear optimal design techniques can be divided 
into two parts: 

1) The selection of the quadratic performance index. 

It can be shown that the selection of the performance 
index can be interpreted as a selection of the closed- 
loop poles of the system in a logical, satisfactory manner. 
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2) Determination of the feedback control law that will 

yield the desired closed-loop pole configuration. 

The optimal design theory approaches a dynamic control problem 
in a reverse sense from that of conventional design procedures. In the op- 
timal approach, the closed-loop poles are selected first and the control 
configuration is then determined. In a conventional design procedure, the 
closed-loop control configuration is first selected, and then the closed-loop 
poles that result from this selection are obtained. 


Summary 

The work presented in this report was motivated by the knowledge that 
use of a quadratic performance index can satisfy many of the control system 
criteria important to the elastic launch vehicle problem. In addition, it was 
felt that the analysis tools developed in Reference 2 would enable a designer 
to systematically specify a performance index that would yield a satisfactory 
control system design. 

This report describes two efficient techniques for optimal launch 
vehicle control system analysis. But, as might be expected, the problems 
associated with the design are primarily ones of synthesis, not analysis. 

The control law must be expressed in terms of measurable quantities, with 
feedback gain magnitudes that do not exceed limits set by previous experience 
in launch vehicle control system designs. The problem is further complicat- 
ed by the fact that the elastic properties of large launch vehicles are not 
known to a high degree of accuracy, requiring a design that is insensitive 
to parameter variations of the elastic characteristics of the vehicle. 

The work described in this report systematically investigates the prob- 
lem areas of elastic launch vehicle control system synthesis outlined above. 

The following section briefly outlines linear optimal control theory in order 
to acquaint the reader with the techniques that are subsequently used in later 
sections of the report. More complete descriptions of the theory are given 
in References 2, 3, and 4. 

The third section describes a method for selecting quadratic perfor- 
mance indices in a way that will yield satisfactory closed-loop dynamics. 

The aim is to provide a suitably fast closed-loop system response and at 
the same time increase the damping ratio of the bending modes. The control 
law is synthesized in terms of measurable quantities. 

In the fourth section, a drift minimum model is used in conjunction with 
optimal control techniques to specify a control system that approximates the 
drift minimum criterion in an error squared sense and within the allowable 
control motions of a single control variable. Direct control law synthesis 
procedures are demonstrated and the optimal control is computed for several 
sets of state variables. A brief investigation was made of the sensitivity of 
the optimal system to variations of the bending mode slopes and shapes and 
preliminary results indicated that the optimal system was remarkably insen- 
sitive to changes in these parameters. The state vector was approximated 
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in a least squared sense by three measurements of the system dynamics. 
Finally, it is shown that higher order structural modes which are not 
accounted for in the optimal control design may be destabilized. There are 
several practical approaches to the solution of this problem. 

These results led to a renewed interest in a study of sensitivity and 
parameter variations, considered in Section 5 of this report. This section 
describes the theory and practice of optimal design of systems subject to 
parameters whose values are known only as random variables described on 
a sample space. A method of optimal compensation is developed in the 
frequency domain and several examples are used to expose the basic analyti- 
cal difficulties which are evolved. The resultant theory is then applied to 
the design of a compensating network required for a flexible booster when the 
value for the slope of the first bending mode is uncertain. 

Section 6 describes a study in which an attempt was made to develop 
a design procedure that limits the feedback gains of one or more state vari- 
ables yet still satisfies a quadratic performance index criterion. A theory 
is developed for this difficult problem and the results show that a direct 
application of the theory requires the solution of a complex set of nonlinear 
algebraic equations. An analog computer program was run to obtain optimal 
solutions experimentally for systems having one or more feedback gains 
arbitrarily set to preselected values. The approach was to start with an 
optimal solution obtained without gain constraints. Then those gains that 
were beyond allowable limits were gradually reduced, while systematically 
adjusting the other gains to keep the performance index as low as possible. 

The cases tried indicated that this procedure is workable and that satisfactory 
closed-loop system designs may be obtained with this approach. 

The report ends with a section on conclusions and recommendations. 
The linear optimal control theory is shown to be applicable to the design 
of launch vehicle control systems. Practical problems of system design 
were considered in the study, such as choice of sensors and their locations, 
sensitivity to parameter uncertainty, and limitations on usable gains. In 
the situations studied, these problems were found to be entirely manageable. 

It was felt that designs of optimal control systems were obtained which met 
practical requirements and which could be mechanized. A systematic ap- 
proach was developed for dealing with uncertainties in the values of system 
parameters. A method was also found for imposing constraints on some of 
the feedback gains and finding optimal values for the system gains within these 
constraints. 

In summary, the study reported here has shown that linear optimal 
control techniques can be used as effective tools in the design of control 
systems for flexible launch vehicles. 
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* 


SECTION 2 

A BRIEF OUTLINE OF 
LINEAR OPTIMAL CONTROL THEORY 


The optimal control problem can be briefly stated as follows: start- 

ing with some initial condition of the plant, the problem is to find a time 
history of control motion that forces the plant to dynamically respond in a 
manner that minimizes an integral function of the control vector and the state 
vector. If the plant is linear and the integral contains quadratic scalar func- 
tions of the state variables and control variables, the control motion can be 
generated as a linear combination of the state variables, and the result is a 
linear feedback control law. This closed-loop system then responds optimally 
to any initial condition of the plant. If the plant is describable by a set of 
constant- coefficient linear differential equations of motion and the integral 
is taken to infinity, the optimal feedback control law is composed of a set 
of constant feedback gains from the state variables of the plant. 


For this problem, the differential equations of motion are assumed 
to have constant coefficients and so can be written as a matrix set of first- 
order equations of the form 

> - F % + Gu Cj = /A: (2-1) 


where v 


- the variables of the differential equations of motion 


^ = the control vector 


{j ~ the output; a transformed set on the state whose 
motions are to be minimized 

F - an n x n matrix of constants describing the coupling 
among state variables in the equations of motion 

G - an n x p matrix of constants describing the effect 
of control inputs on the equations of motion. 


A control motion a & 
performance index 


is to be found that minimizes the quadratic 
oo 


2 V 7 =• J (if' Qy +■ u' f2 u.)dt 


( 2 - 2 ) 


where Q - an r x r positive definite symmetric matrix whose elements 
weight the contributions of each output in the integral 


13 = a p x p positive definite symmetric matrix whose elements 

weight the contributions of each control motion in the 
integral. 


There are several ways to obtain the optimal control motion u Q that 
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A 


minimizes the performance index. In the time domain, the calculus of vari- 
ations may be used to obtain a solution (see References 11 and 2). If this 
approach is taken, the Euler-Lagrange partial differential equations of mo- 
tion must be satisfied. 



d 

(M\ 

= 0 


dt 

[St) 

dx 

±i 

! d%\ 

= 0 

da 

dt ( 

9&J 



where <X is the Lagrangian 

X - { (y'Qy +u'ieu)* + r* + Gu) 

and % is an undetermined vector called the adjoint state vector or the 
costate. 


The Euler-Lagrange equations, together with the original equa- 
tions of motion are 


% - - & u 

= 0 

(a) 


ftu 0 + G ' % 

= 0 

(b) 

(2-4) 

hi ' Q 1-1% + % *■ F % 

= 0 

(c) 



From Equation 2-4b the optimal control law is u 0 = -P’ 1 G*%- 

If the substitution % - P X, % - Pi is made, Equations 2-4 can be 
reduced to the single quadratic Riccati equation in P which was described 
by Kalman (Reference 3) and others. 

0 * VF + F'P-VGi1Z-*G 9 P + H 9 Qil (2-5) 

The feedback gains are determined from the positive definite symmetric 
solution of the Riccati equation. 

Now consider the Laplace transform of Equations 2-4: 


~ Is-F 

-G 0 


^( 5 ) 


%(o) 

0 

ft G’ 



- 

0 

-H'QH 

0 -Is-F' 

1 

%(s) 


-%(o) 


The determinant of Equation 2-6, when set to zero, yields an expres- 
sion for the closed-loop roots of the optimal system with state variable v* and 
"adjoint" system with adjoint variable % . The roots of this characteristic 
equation can be determined from the root square locus expression 

| j. e- , a'[-is-rT , //w[is-rY / a I = o (2-7) 
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which can be derived from the determinant of Equation 2-6 (References 2 
and 5), or from a frequency domain formulation of the problem (Refer- 
ences 2, 9). H[ls- is a matrix of transfer functions ^/ Uj ($) relating 

the outputs of the system to the inputs and 'H* is the transpose of 

the transfer function matrix with s replaced by - s. 


The linear optimal control problem can also be solved entirely in 
the frequency domain. Parseval's theorem is applied to the performance 
index. Equation 2-2; •<*> 


where 


2V -- 


/ 

2/rj 


(y^Qy + u„12u)ds 


(2-8) 


y - y (s) = /■/ \l s - F~\ G u. + H 5 - /- 1 x ( o ) 


= YJu, + 3x(o) 

y* = y'(- s (o) S„ 

% - *'(s) 

W* = 


It is desired to take a variation to find the control that minimizes 
2-8. To do this, substitute for (j and y* in Equation 2-8 and let 

u, = d a + 


It will be found that the requirement for minimization (References 
2 and 9) is that j o o 

% ^ vJu, a + W* <3 (o) i'Bu.^ds = 0 


X. 


ITT; 




or 


["£? +- IV# Q Vs/] u, a +■ W# Q ’B'X (o) - ( s j 


(2-9) 


where ^ ( 5 ) must be analytic in the left-half plane. * Equation 2-9 is a 
matrix equation of the Wiener-Hopf type that must be satisfied to obtain the 
optimal control vector. The solution to Equation 2-9 is given by 


( 2 - 10 ) 


where 

and Y,-'[W„08*M] 
sion into the sum 


Y Y* - TlS + W„(3W] 

has been decomposed by a partial fraction expan- 


c(o)] = [y.'V]^ [v,"Y. 


( 2 - 11 ) 


* Appendix B describes this requirement when the open-loop plant is non- 
minimum phase. 
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It can be seen that 


1 13 + W* Q W) = 0 (2-12) 

yields the same expression for the root square locus as Equation 2-7. 

A complete method of analysis and synthesis is therefore obtained. 

The selection of a performance index that yields a satisfactory closed-loop 
response is possible through the use of the root square locus expression. 

The feedback control law can be obtained by solving the Riccati equation 
or by satisfying the matrix Wiener-Hopf equation. References 2 and 5 de- 
scribe the solution through the use of the Riccati equation, while References 
2 and 9 describe the proper technique for solving linear optimal control 
problems using the matrix Wiener-Hopf equation. 

When the plant has only one control input, the synthesis problem can 
be solved in a direct, straightforward manner. The plant is described by 
the set of first-order linear equations of motion 

i, = F% + Gcl 

and the optimal feedback control law is of the form 

U 0 (2-13) 

The equations of motion of the optimal regulator are given by 

% = (F- (2-14) 

whose closed-loop characteristic polynomial is 

Ms) «|ls-F +Gkt\ (2-15) 

The closed-loop character! Stic polynomial can be reconstructed for any 
values of weighting in the performance index from the root square locus 
plot. Equating powers of s of the characteristic polynomial obtained from the 
root square locus plot to the scalar expression of Equation 2-15 yields a 
linear algebraic set of equations in the k 4 - feedback gains. This set of linear 
algebraic equations can be simply solved by machine in a straightforward 
manner. Examples of optimal analysis using the root square locus plot and 
synthesis using the method described above are given in the next two sec- 
tions. In addition to the examples in this report, Reference 2 contains num- 
erous examples designed to acquaint the engineer with the required mathematics 
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SECTION 3 

SELECTION OF THE PERFORMANCE INDEX 


Introduction 


The selection of the parameters of the performance index is not 
automatic. The variable or variables that are included in the performance 
index will have an optimal response that tends towards the response of a 
Butterworth filter as these variables are weighted more and more heavily 
with respect to the control. The implications of this statement are signifi- 
cant. It means, for instance, that if the performance index contained only 
a bending mode variable, say , then would tend to respond as a 

second-order system with a damping ratio of approximately % - 0. 7. The 

responses of the other variables, however, are almost certain to be unaccept- 
able. The approach, then, will be to include a measure of both a rigid-body 
variable and bending mode variables in the performance index. This approach 
will permit a rapid and stable response of the rigid-body contribution to the 
dynamics and acceptable damping of the bending mode contribution to the 
closed-loop response. 


One Bending Mode 


Although it is clear that a measure of both the rigid-body mode and 
the body bending modes will be included in the performance index descrip- 
tion, the mix of the two variables should be investigated. Thus the attitude 
angle and first bending mode variable may be included in the performance 
index in the following ways: 


00 


2V - 


'min 

fit. 




(3-1) 


or 


00 




(3-2) 


The analysis of the first performance index, Equation 3-1, is a rela- 
tively simple matter, involving only a single root square locus plot with 
the parameter ql f . A preliminary analysis of a conventional root locus 
plot of the function . 

. 08 , , ft . . „ 

< 3 - 3 > 


will yield the zeros of the function that will appear in the root square locus 
plot. In general, the roots of Equation 3-3 are chosen to obtain an adequate 
separation between the poles and zeros of the function of Equation 3-3. This 
will guarantee that adequate damping of the bending mode will be achievable 
for some value of the performance index parameter 


The analysis of the performance index of Equation 3-2 involves two 
root square locus plots instead of one conventional root locus plot and one 
root square locus. As before, the first locus serves to define the zeros 
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for the second root square locus plot. The root square locus expression 
for the performance index of Equation 3-2 is given by 

*» ' ' 


1 + — (*) (s) 

r /3 C (3 C 

where the term in brackets 


/ + 


^ z fic 


1± 

(3c. 


Lh)^(-s) 


} ' “U 


(-sj 


» 0 


(3-4) 


1 + 




g f5) a ( - s) 


= 0 


(3-5) 


defines the zeros, or terminal points for the locus of the closed-loop poles 
of the optimal system. Because of the quadratic nature of Equation 3-5, 
the locus of the zeros of Equation 3-5 will generally have a higher damping 
ratio for any value of the parameter q 2 j </ t than will the zeros defined by 
Equation 3-3 for any value of the parameter k. 

To illustrate this process of the selection of the terminal points of 
the root square locus process, it was decided to use the performance index 
of Equation 3-1 for the analysis of an optimal system that contains only one 
bending mode and the performance index of Equation 3-2 for the two bending 
mode analysis and synthesis problem of this section of the report. 


Single Bending Mode Analysis 

The elastic booster equations of motion including a single bending 
mode are (see Appendix A): 


0 e - .0733 cc + 0.45/3= 0 

+-.D405 6t + . OW67oc, + .OZH (3 = 0 
- 545 3 oc + + . 023 Z ^ * 5. 37 - (5-83/3= 0 

(3 + 1 7.9/3 = 17.9 /S^ 


(3-6) 


or, in matrix form 


— 1 
eo 

l 


0 

1.0 

0 

0 

0 

0 


I 

1 


0 



0 

0 

.0733 

0 

0 

-.45 




0 

« 

OC 

— 

-. 04o5 

1.0 

-.0107 

0 

0 

-.0211 


CC 

+ 

0 



0 

0 

0 

0 

1.0 

0 


v, 


0 



0 

0 

5.453 

-5.37 

-. 0232 

15.93 


i 


0 



0 

0 

0 

0 

0 

-17.9 




17 . 9 


fit 

(3-7) 
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The design will involve the analysis of the performance index 

00 

m x /[?(?, ' * r p c *]<* 


(3-8) 


It has been shown in Section 2 that the closed-loop roots of the optimal 
system and adjoint are given by the root square locus expression 


where, for this example 

9-M 


= o 


(3-9) 


//- [H 


o 0 


P [Is-F]' £ = 


/ o o 
n 


] 


T2' f = 


r 


, 

TP 








Pc 


A preliminary root locus plot of ( s) + £ — ^ (s) = 0 was performed to 

pc* pc 

select a suitable value for k. Since this locus defines the zeros of the root 
square locus plot, the selection of k was made on the basis that the roots of the 
closed-loop optimal system would be such to yield a sufficiently rapid booster 
response and increase the damping ratio of the bending mode. A value of 
k = 30 produced the desired result. This locus is not shown because its con- 
struction is straightforward and relatively simple. 

Substituting the transfer functions from Appendix A, 

n t 0* 

V(s) = — (s) + 30— (s) 

Pc PC 


<*>« 



Substituting P(s) and P(-s) into the root square locus expression (Equation - 
3-9) yields 


'4 


3027 ft 


.0H5/ 


1 ± 


2.29, 


ft 


2.24 


It 


.0418 


It 


.242 


ft 


294 


ft 


I7.9J 


J 2 (.005) _ 

ft — 5 s + 

1 2.317 


—Y 

'■■VlJ 


.(3-11) 
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One half of the solutions to Equation 3-11 are plotted in Figure 3.1. 
The stable, optimal solutions are shown and the unstable, adjoint system so- 
lutions are omitted for simplicity and clarity. The parameter of the locus is 
q/r, the ratio of the weighting of the state variables and the control variables 
in the performance index. The locus then shows the closed-loop poles of the 
optimal system as a function of the ratio of the weighting of the two parts of 
the performance index. For purposes of illustration, a value of q/r = 0.01 
was chosen as a candidate for a weighting of the performance index factors 
yielding satisfactory closed-loop optimal system dynamic characteristics. 

From the root square locus plot, a value of q/r = 0.01 yields the 
closed-loop characteristic polynomial 

A(s) = (s*f7.9)(s + .0//7)[_S z +2(.8f)(t./3)s + Z(lz)(ZAZ)s + (2.nf\ (3-12) 


Equation 3-12 shows that the damping ratio of the first bending mode has 
been increased from an open-loop value of ? = . 005 to ? = . 12 in the closed- 
loop design. At the same time, the rigid-body contribution to the response 
contains a predominant pair of pole s at cd n = 1 . 1 3 rad / sec, % =0.81. The at- 
titude response of the vehicle will therefore be smooth, with little overshoot. 

The optimal control law that yields the closed-loop characteristic equa- 
tion 3-11 and satisfies the performance index 

oo 

(3 ' 13) 

0 

can be obtained by solving the Wiener-Hopf equation, 2-9, or the matrix Ric- 
cati equation, 2-5, or by a direct method that is described in the next section. 
For this problem, the feedback control law is found to be 

/% = 3.00 +4.980 e + +A5aL-.0827f t - -0Z3lh -0.135/3 (3-14) 

The magnitudes of the feedback gains are generally within the guidelines 
outlined by the Marshall Space Flight Center of NASA, and can be mechanized 
in terms of measurable quantities. However, it is clear that the quantities con- 
tained in the optimal control law of Equation 3-14 do not generally represent 
measurable quantities. It would be more realistic to express the feedback con- 
trol law in terms of the outputs of sensors located on the vehicle. Because 
this particular problem, which includes only one bending mode, has been pre- 
sented for illustrative purposes only, the feedback control law mechanization 
aspects are not presented. The remainder of Section 3 and all of Section 4, 
which includes a description of the vehicle containing two bending modes, does 
consider control law synthesis problems in detail. 

Selection of Performance Index for Two Bending Modes 

The addition of a second bending mode to the description of the launch 
vehicle adds complexity to the problem because now a measure of the second 
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Figure 3. 1 Locus of Poles of the Optimal System 
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mode dynamics must be included in the performance index. The equations of 
motion for this problem are (see Appendix A) 


1 «« 

1 


* # 

0 

ne 


aC 

- 

i 


% 




h 


i 



0 10 0 

0 0 .0737 0 

.0405 I -.0107 0 

0 0 0 0 


0 0 

0 0 


0 

-.45 


0 0 -.0211 oO 


0 0 0 0 / 0 

0 0 545 -5.37 -.056 3 0 

0 0 0 0 0 0 

0 0 2.36 0 0 -31.1 


0 15.6 3 q 


fit 

( 3 - 15 ) 


0 -31.8 -.05 64 22.77 

0 0 0 -1 7.9 


A performance index is chosen of the form 

OO 


zv - IV b, 


(3-16) 


requiring the selection of three weighting parameters of the performance 
index. The analysis of this problem can be carried out using the root square 
locus expression 




( 3 - 17 ) 


where for this problem 

| 1 0 0 o o o o o 


A/ = 0 0 O 1 0 O 0 O 

0 0 O 0 O I o o 


% 0 0 

Q - 0 ° 

0 0 9z 


= F 


( 3 - 18 ) 


The expr es si on H [i S - F ] 1 G is a column matrix of transfer func- 
tions of those variables appearing in the performance index. 








liSil ' EMil 







IH-2089-F-1 


Substituting into Equation 3-17, 



— (-S 


n t . , 

— (-s) 


o 

1 

0 

Q, 

0 

0 

r 

fie 


ft J 


0 

'3 

0 



or 


n 4 0 *, , $*, , 9 2 n 

0 — / + — fs) — (-sj + 


r fic 

Equation 3-21 can be written 


0 = / 4 -E -A(S) A(-Sj < /* 


(*) 

r fie pc 


k^ 

A<s) I (3-20) 

> 


*> *L(,)*L(. S ) (3-21) 

fie fic 


r fit (3c 


fi _ * fit (gj fie C ' 5) 




Qz ^2 / \ 

lf (s ^ -tt ^ 5 

pc /3c. 


(3-22) 


Equation 3-21 shows that the closed-loop poles of the optimal system 
and its adjoint are a function of three parameters q t (r , , and <f. 3 jv . 

Equation 3-22 shows that three separate root square loci are involved in the 
analysis of the closed-loop system and the selection of the values of o If. 
%/r > and qjf. ^ 


Defining 


n ’ („ k (.,) , KL 2 


fi. 


n 

~fic 


ftc 

h 

Pc 


1 (s) (-S) = 


— ( 3 ) -^(- 5 ) » — 

fie fic 0 D & 

then the first of the three root square loci can be written 

0-1 + — 


"ft 

DP 

DP 

^0 ^0 


(3-23) 


V 


/V, Ay, 


(3-24) 


Substituting the numerical values for the booster transfer functions, Equa- 
tion 3-24 becomes 


0 -, + 2* <*">’ 


.0408 


.454/1 .499 


5.639 


5.639/ J 


9* it 


.0412 


.3/9 


s \/ s 

H 


.369 


J 2 Coos) 

1 ± 5 + 

2.317 


2.317/ 


(3-25) 
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One half of the root square locus of Equation 3-25 is plotted in Figure 
3.2, The stable, optimal system part of the locus is shown, and the unstable, 
adjoint part is omitted. The locus in itself is difficult to interpret precisely, 
but a few general comments can be made. This locus and the one to follow 
will define the zeros, or terminal locations for the closed-loop poles of the 
optimal system. Therefore, in general it can be said that the aim of the 
locus of Figure 3.2 will be to select roots that have an adequate damping 
ratio. This will guarantee that both the first and second bending modes will 
have adequate damping ratios in the closed-loop design. The locus does indi- 
cate that if good first mode damping is to be achieved, should not have 

too small a value, while too large a value for will result in an insuffi- 

cient second mode damping ratio. A selection of 4 ? / q, z in the range 
. 2< “?3 / <2 will likely serve the purpose of obtaining good damping ratios 

for both the first and second modes. 


The value ^ 3/^1 = 0* 5 was chosen as a logical value, yielding a max- 

imum damping ratio for the roots of Figure 3.2. With a value for q. 3 / q 2 
chosen, the second root square locus expression to be formulated is 


o =1, 


s 


_2t ' - 


% 


k? 2 


(Roots from Figure 3.2 with - .05) 

_ 


(3-26) 


Substituting the appropriate dynamic quantities from the transfer func- 
tions and from the first root square locus expression, the second root square 
locus expression becomes 


0 


1* 9.SZ 


% 

% 



2 639 ; 
3.92 




2 (. 005 ) / s \2 

2.3/7 5 + (2.3/7/ 


2 693 ) 
.47 




Z(.Q05) /_£_) 

S .& 39 S \5.6> 39/ 


2 


(3-27) 


As in Figure 3. 2, one half of the root square locus expression of Equa- 
tion 3-27 is plotted in Figure 3.3. This locus defines the end points or zeros 
of the root square locus for the closed-loop poles of the optimal system. An 
inspection of the plot shows that a value of q z = 4 allows for good damping 
of both the first and second bending modes. 


The third root square locus expression can then be formulated as 


0 = 




(Roots from Figure 3.3 with 


4) 


DD 


(3-28) 


where_ D is the open-loop characteristic polynomial of the elastic vehicle 
and D is the open-loop characteristic polynomial with s replaced by -s. 


Substituting the proper expressions from the second root square locus 
plot and from the transfer functions yields 
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igure 3. 2 Root Square Locus Plot of Equation 3-25 




igure 3. 3 Root Square Locus Plot of Equation 3-27 
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The root square locus plot of the poles of the optimum system is shown 
in Figure 3.4. As in the previous plots, only the stable, optimal half of the 
locus is shown and the adjoint system is omitted. The locus shows that the 
objective of the design procedure has been achieved. The open-loop poles of 
the booster originating from the rigid-body dynamics become stable, well 
damped and increase in natural frequency as the parameter of the locus 
is increased. In addition, the damping ratios of the poles originally associated 
with the bending modes increase significantly. The net result is a logical de- 
sign procedure that yields a desirable and satisfactory pattern of closed-loop 
poles. 

Optimal Control Law 


It is a relatively simple matter, once the root square locus plot has 
been obtained, to calculate the optimal control law that will yield the closed- 
loop polynomial selected from the root locus plot. The original open-loop 
equations of motion are written 

% = F % + Gti 

and it is desired to find the optimal feedback control law 

u 0 - - 

so that the closed-loop optimal system becomes 

£ - (F-GtS)x 

The closed-loop characteristic equation is given by 

I Is- F + Gh 1 1 = 0 (3-30) 


Equating the determinant of Equation 3-30 to the closed-loop polynomial ob- 
tained from the root square locus plot yields a set of linear algebraic equa- 
tions from which the feedback gains can be easily computed. 

As an example, a value of Cf tf / f = 50 was chosen to yield a closed- 
loop system whose speed of response would be satisfactory. From the root 
square locus plot of the poles of the closed-loop optimal system (Figure 3.4), 
the closed-loop polynomial is given by 


h(s) = ($+.08Z5)(s+t7.9f)[s Z 4-2(.76)(!.Z8)s-+(/'28)*\[sh2(.t7)(2.4)s + (2-4) 2 ] \s z *2(.06)(f.66)s* (&(>5) Z ~\ ^ 3 1 ^ 
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For this problem, the determinant A (s) — + or 


A(s) = 


S 

0 

.0405 

0 


- 1.0 


-.0733 


1.0 s+ .0107 


0 


0 


-6A53 


0 


-2.36 


5137 


0 


0 

0 

- 1.0 

S+.0232 


0 

0 

0 


0 

0 


0 

O 

0 

0 

0 

- 1.0 


0 

.45 

.0211 

0 

-15.23 


0 


31.8 $+.0564 -22.77 


17. 9jfe f /7.9 17.9t 3 17.9 17.9H S /7-9^ 6 !7.9t 7 5+17.9+ 17.9t a 

(3-32) 

Equating powers of s of Equations 3-31 and 3-32 and solving for K yields 

\5 = [~ 9. 104 , -5.993, 4-407, .1020, .05163, .05513, .02264, . /86o\ (3_33) 

and the optimum control law becomes 

(3-34) 

/% ^ +9.104 fy+S'W3(p e --4'.4t>7cC--1020y f -'03f63)2 f '.<)3Si3y z -.02Z6% z -. 186/3 


The control law of Equation 3-34 is not realizable because the variables 
of Equation 3-34 are not directly measurable. It is necessary to define mea- 
surable quantities and express the control law in terms of these quantities. 

For purposes of illustration, the following set of state variables was chosen 




9 ^ 3 



(3-35) 


where 0~, , 0^ . , and represent the quantities measured by a 

1 & s V i 

position gyro, a rate gyro and an accelerometer at station ^ on the launch 
vehicle body. Choosing the body stations X f = 41.5 meters, % % = 86.0 

meters, and = 122.5 meters, the transformation =? A lo becomes, 

for this example * 
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.060 
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*1 
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0 

0 
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.01624 

-35.38 
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3.550 


I 

% 

K 


0 

10 

0 

0 

.1230 

0 

-.143 
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% 

% 


0 

0 

23.3 7 

-14.32 

-.0473 

-8.526 

-.01861 

14.287 


i 

j _ 


0 

0 

0 

0 

0 

0 

0 

1.0 




and the control law becomes 


A = 9.m0„ +4.2260„ -.4202A* +t 3)7 <ZL -.1770.* + .4-53 <$#-. 0969 a, -.1493/3 

1 ° 7 7 AJfc ^ /Jfg 5 ^ 


(3-37) 


sensors located sensors located sensors located 

at station 41.5m at station 86.0m at station 122.5m 


The control law of Equation 3-37 represents a solution to the two 
bending mode problem. By locating three instrument packages at three dif- 
ferent locations along the body of the launch vehicle, every state variable of 
the system is measurable. The feedback gains are low enough that the con- 
trol law may be mechanized without difficulty, resulting in a design that may 
be practical for a launch vehicle whose significant elastic properties can be 
described by two bending modes. 

The transient responses of the optimal system, defined by the control 
law of either Equation 3-34 or 3-37, are shown in Figures 3. 5 through 3. 7. 
The characteristic polynomial of the closed-loop system is given by Equation 
3-31, which indicates that the first mode damping has been increased from 
X = . 005 open loop to approximately X = . 17 for the closed-loop design, 
while the second mode damping ratio has been increased by a factor of 10, 
from X = .005 to approximately X> = .05. The transient responses to a unit 
step command input clearly show the increases in the first and second mode 
damping ratios. In addition, the rigid-body contribution to the closed-loop 
response, dominated by a pair of poles located at c£> n = 1. 28 rad/sec and 
X = 0.76, contribute to the smoothness of the vehicle’s response. The con- 
trol motions of Figure 3. 7 show that the gimbal deflections are reasonable 
and well within the frequency capability of the actuator. 
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Figure 3.5 Normal Acceleration Response at 
Three Different Body Stations 
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SECTION 4 

MODEL ORIENTED DESIGN STUDY 


Introduction 


This section describes the use of a model in the linear optimal con- 
trol design procedure. The model defines a set of equations that describe 
the motions of an ideal launch vehicle. Optimal control is then used to se- 
lect a feedback control system configuration that minimizes the differences 
between the motions of the launch vehicle and the model in an integral error 
squared sense and control squared sense. Minimizing the motions in this 
manner guarantees a stable solution and also yields an error response that 
is smooth and well behaved, resulting in control motions that do not require 
violent deflections if the model is smooth. A well behaved model is desirable 
for the launch vehicle, so it is expected that the control motion will be 
smooth for this application. 


The following theoretical development and examples show that the use 
of a model is a logical design aid, resulting in a more direct and simple anal- 
ysis procedure than the procedure presented in Section 3. The optimal aspect 
of the design tends to satisfy some of the design objectives, such as stability 
and smoothness of response, while the model aspect satisfies other require- 
ments such as drift minimum characteristics and speed of response of the 
closed-loop system. 


A model of a desirably responding system is formulated as (Ref. 3) 




(4-1) 


where L defines the matrix of constants describing the coupling among vari- 
ables in the equations of motion of the model. 


A control motion u 0 is to be found that minimizes the quadratic 
performance index 


ZV = f Jy-Ly)'Q(y-Ly) + u. , '&u\a / £ 


(4-2) 


where 


Q = an rxr positive definite symmetric matrix whose elements 
weight the contribution of each error in the performance index. 

V, - a pxp positive definite symmetric matrix whose elements weight 
the contribution of each control motion in the integral. 

The error portion of the performance index, the {y-Ly) term, is a 
vector that vanishes when the closed-loop optimal system behaves exactly 
like the model. In general, the feedback cannot exactly match the booster 
to the model for several reasons. In general, the L and F matrices are 
of different dimension, and a limited number of controllers is available for 
feedback purposes. 
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Although it will be impossible to exactly match the booster and the 
model, the error can be made as small as practical by weighting the error 
portion of the performance index heavily with respect to the control. It will 
be shown that approximating the model response through the use of a qua- 
dratic performance index results in not only an optimal, but an acceptable 
feedback control system design. 

The problem of minimizing Equation 4-2 subject to the differential 
constraint of the equations of motion is again a straightforward problem in 
the calculus of variations. The Lagrangian of this problem can be defined 
as: 

cC “ \jy -^y)'Q (q-Ly) + u' + %'{- y ± Fy. + Gu'j 

-{ \y' IJ'QUY'- vF 'QLHp-v'F'L’QF* (4 _ 3) 

-f yJL' hl'QLH% + u + %'(- y ■* Fy. + Gu) 

where % is an n x 1 vector called the adjoint state vector or the costate. 


The minimization of the quadratic performance index requires that 
the Euler -Lagrange equations be satisfied 


dX d_ 1 9cC_ ) _ 

2 y ~ dt [Gy,/ 


= 0 


2X d IdcC 


= o 


(4-4) 


are 


du. dt ( 9d 

Using the Lagrangian of Equation 4-3, the Euler-Lagrange equations 

% + F%-H'QM + U , QLU* -F'L'Qdy +F'L'QLFy = 0 

(a) 

(4-5) 

(b) 


12u,„ + G'%-0 


The optimal control law from Equation 4-5b is u 6 - -F' 1 G r % and 
the optimal closed-loop system becomes 


y, - F% + G1Z' 1 G > % = 0 


(4-6) 


The exact solution for % as a function of % has been obtained 

by Kalman (Ref. 3) and Tyler (Ref 6). For the case of the single control 
input the feedback control law can be more quickly and easily obtained by 
developing a root square locus expression, which, when plotted, spectral 
factors the poles of the closed-loop optimal and adjoint systems. The closed- 
loop characteristic polynomial of the optimal system can then be recon- 
structed from the root square locus plot, and the feedback gains are then 
easily obtained. The root square locus expression is developed below. 


In Laplace transform form, Equations 4-5a and 4-6 become 


[Is- A] 

GK-’G’ 



%(o) 

-//'[-Ts-2/]Q[ls-L]^ 

tii-r] 

%(s) 


*(o) + v(oj] 

_ -d'\QL +L'<d[ dv(o) 
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The determinant of the left hand side of Equation 4-7, when set to zero, de- 
fines the characteristic equation of the optimal system and its adjoint. It can 
be shown (Ref. 2) that this determinant may be written 


I \ls-p] GG-'G' 

-//'[-Is- [is - F'] 


* je|ls-F||-Is-r| I^- , fi , [-Is-F']'V'[js-/.'].(?[/s-Z.]w[ls-F]' 



(4-8) 


if [is-r] an d[-Js-r'] are square, non-singular matrices. Rut the determi- 
nants /ls-F| and j-Is~ F r l are the characteristic polynomials of the open- 
loop system and its adjoint, which vanish only at their root locations. The 
closed-loop roots of the optimal system and adjoint are therefore obtained 
from the expression 


|l4 IS-tG'l-Ts-F'yU'l-Is-L^Q [ls-L]hl[ls-F]’ f G | = 0 ( 4 -9) 

The term A/ [js - F~\ is a matrix of transfer functions (s) relating 

the outputs of the system to the inputs and Js-F'] /-/' is the trans- 
pose of the transfer function matrix with s replaced by -s. The locus of 
roots obtained from Equation 4-9 defines a root square locus expression 
(Ref. 2 and 5). Using this expression one can investigate the closed-loop 
characteristic polynomial of the optimal system as a function of the weight- 
ing parameters of the performance index of Equation 4-2. Because the root 
square locus is a p x p expression, while the original determinant of Equa- 
tion 4-8 is of dimension 2n x 2n, a considerable simplification has been made 
for purposes of analysis. In particular, because the booster has a swingle 
controller (the gimballed rocket engines), the root square locus reduces to a 
scalar expression. 


Basic Limitations of Modeling 


The use of a ’’mathematical model” to define a control system cri- 
terion has limitations and clearly cannot be used without some conception 
of the implications of its use. In general, it is mathematically and physi- 
cally impossible to exactly match a system to a model. Generally speaking, 
there are two reasons for this: 

1. The mathematical model is usually of an order or dimension that 
differs from that of the actual system and often differs from the 
mathematical description of the actual system. 

2. In general, the number of controllers does not equal the number 
of degrees of freedom of motion of the system. It is basic 

that the number of controllers must equal the number of degrees 
of freedom the plant possesses to have complete control of the 
dynamics of the plant. 

To illustrate the two principles stated above, consider the following 
simple examples: 

1. Assume a second-order, single controller plant is describable by the 
first-order set of equations 

'V' - F % 4- <3u* 
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— -i 


r- 



*,t 
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> 
1 

1 

> 

$22 


% 2 


1 

1 


Assume that it is desirable to use feedback control of the form 

a = -1^ X 

in an attempt to match the plant to a model described by the matrix 


(4-10) 


(4-11) 


L 



12 


I 


Z1 


1. 


22 


(4-12) 


Substituting Equation 4-11 into 4-10 yields, in the closed loop, 



r ~ 




— 


— 



> 

) 

r~ 






$12 



9i ^ z 

1 

1 

) 

*1 

(4-13 


*2 



$21 






f 

% 2 


or 

L J 


\ 


-J 


— 

— 1 

I 


- - 



% = [F-Gtf] % 

It is clear, by an examination of 4-13 that only one row of L, either 
{_^zt ^21 1 or [ A/ &iz ] can k e matched by the two feedback gains 

and k^. 


The best that one can do is to specify L and F as companion matrices, 

i. e. , 



— 

— 


— 

— 


0 

1 


0 

i 

L - 

z 

~ 2 

F = 




Then feedback control can be used to match the natural frequency 
and damping ratio of the plant and the model. This approach has been incor- 
porated into the design philosophy of the flexible launch vehicle. 
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2. Assume that the plant is describable by the third-order set of equations 


’ V 


*t 1 

*tz 





V 

^2 

= 

*2, 

^ 22 



*2 

+ 


_ V 


_^3f 

^32 



_*jJ 


_ ^3 _ 


and the desired "model" is described by the same 


2x2 matrix 


L 


p 11 


JL 


21 


l 


12 


JL. 


22 


It can then be deduced, using the principles from example 1 that it 
is possible to match only one row of any 2x2 submatrix of F to L. If F and 
L are specified as companion matrices, then the two poles of the model may 
be exactly obtained and the third pole "arbitrarily" selected. 


The approach to flexible launch vehicle control system design in- 
corporates both of the principles described above with one additional, and 
perhaps important, factor. Referring to example 2, one may ask the follow- 
ing: 

Instead of exactly matching two of the model poles and "arbitrarily" 
selecting the third pole, what distribution of the three (closed-loop) 
system poles approximates the two model poles in a minimum 
quadratic "error" squared and "control" squared sense? 

The reason for specifying the closed-loop poles in this sense is logi- 
cal. It eliminates the arbitrariness of the selection of the third pole. For 
if the third pole were too "large", the control motions would be rapid and 
high feedback gains would probably be required. If the third pole were chosen 
too "small", it would contribute significantly to the closed-loop system re- 
sponse. Therefore, the quadratic type of design used in this report is a com- 
promise design. 

If, in the flexible launch vehicle design, the exact drift minimum prop- 
erty must be obtained, it would be a simple matter to precalculate the feed- 
back control required to yield a pole at the origin, then select the remaining 
poles in a minimum quadratic "error" squared plus "control" squared sense. 
The penalty incurred is that greater or sharper control motions would be 
required than for the case when all the plant poles were chosen to satisfy 
the quadratic performance index criterion. Perhaps the penalty is significant; 
more likely the differences would be negligible. 


Elastic Booster Application 


Because the booster has a single control element, feedback control 
can alter only the characteristic polynomial of the vehicle and not the zeros 
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of the transfer functions (unless other dynamic elements are deliberately added 
to the system). Therefore, the optimal control problem may be formulated 
to express the error in terms of a difference between the characteristic poly- 
nomial of the elastic booster and the characteristic polynomial of the model. 

To do this, assume that the F matrix associated with the booster is of dimen- 
sion n x n and L is dimension Ji * £ . Define H of dimension 1 x n 


N = 


K o 

0 k 


n 


0 0 


o'h u 0 


(4-14) 


where k ff = k Z2 = - . . = k tg = / 

The matrices F, L and G are written 


F = 


o 

o 

o 


0 1 0 ... 0 


0 

0 0 1 


0 0 1 


0 

• 

• 


• • 


• 

. 

L * 

• » 

Q = 

• 



• » 


• 

1 0 


/ 0 


• 

0 0 / 


0 1 


0 

•*>» ' • • ’ h n -1 - b » 


'df - d% • • • ~ d £ i -d £ 


1 


(4-15) 


where 


„n / n-> i rt-2 

S + b„ S + b„. f S 4- 


$*+ d £ s 1 '* +dj_ t s £ ' 2 + 


■ b 2 s 4- b 1 = I Is - F I 


. . + d z s 4-d f = 1 1 s-L 


(a) 

(b) 


(4-16) 


The Q matrix is defined as Q - ql , where I is an Ji x identity 
matrix and 12 - T , a scalar. 

When the quantities F, G, H, L, Q, and R are substituted into the 
root square locus expression of Equation 4-9 one obtains 


1 _ q I Is-Ll (-Is-/// 


(4-17) 


Ils-FH-Is-F'l 
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The numerator of Equation 4-17 contains the characteristic polynomial 
of the model and adjoint while the denominator of Equation 4-17 contains the 
characteristic polynomial of the booster and its adjoint. The locus of roots 
of Equation 4-17 will therefore originate at the launch vehicle open-loop 
poles and terminate at the model singularities. The parameter of the locus 
is proportional to q/r, the ratio of the weighting of the error to the control 
in the performance index. The net result of the use of Equation 4-17 is that 
the performance index has been formulated to minimize the difference be- 
tween the characteristic polynomial of the booster and the model in the inte- 
i gral error and control squared sense. 

Equations of Motion of the Elastic Booster 

I ~ r_ " ~ — 

The equations of motion are taken from NASA f, Model Vehicle No. 2 n 
(Appendix A) which is representative of the post-Saturn vehicles under con- 
sideration by NASA. Assuming complete vehicle symmetry, the equations 
of motion in the x- ^ plane at t = 78 sec after launch (max q) are given below 
with certain simplifications. Two normal bending modes are included and 
, the engine actuator is assumed to be described by a first-order differential 

! equation 


0s -0.0733CC + 0.45/3= 0 
-&Z + 0. 0405 0 e + 01067 oc + 0. . 211/3 - 0 

j8+ 17,9/3 * 17.9/8 C = 0 

- 5.4530c + 0.02317 + 5.37 -1563/3 = 0 

- 2 . 360 & + $ f 0.05642 j z + 31.80 y/ z - Zl. 77/3 - 0 


pitching equation 
heaving equation 

actuator dynamics (4-18) 

I s ! normal bending mode equation 
2 n d normal bending mode equation 


where 


0# = rigid-body attitude angle 

= angle of attack 
/3 = control deflection angle 

= first normal bending mode variable 
% 2 = second normal bending mode variable 


The F and G matrices are written in the first-order form of Equation 

4-15 


% ~ F % +- Gu y - A% 

and the output matrix A is defined to be an 8x8 matrix of constants 


31 



IH-2089-F -1 


A « 


a 

a. 


a„ 


a i2 & ti 


'll 


a 


in 


a 


in 


(ti-i)n 


a n(n-l) a »» 


(4-19) 


Then 


1l 

LI 


0.- S 
in 


n -1 Yt-Z 

a in-1 s * 


a., 


- n , YI-1 , L 
$ + bfl 5 


5 ri ' Z -h 


bi 


n 


is the transfer function 


relating the </ t - ^ variable of the equations of motion to the control input. 

In terms of the numbers of Equation 4-18, the first-order equations become: 


if 


0 

1 

0 

0 

0 

0 

0 

0 



r 

0 



0 

0 

1 

0 

0 

0 

0 

0 




0 

is 


0 

0 

0 

1 

0 

0 

0 

0 




0 

if 


0 

0 

0 

0 

1 

0 

0 

0 




0 



0 

0 

0 

0 

0 

1 

0 

0 


% 

+ 

0 

H 


0 

0 

0 

0 

0 

0 

i 

0 


^6 


0 

K 


0 

0 

0 

0 

0 

0 

0 

1 


*7 


0 

*8 


-9.074 

223.5 

-20.71 

-3010 

-193.7 

-(,(,5.5 

-38.11 

-/799 




1 


fit 

(4-20) 
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Cj = AV 




-19.4 

-1376 

-12 A 

-29 9.4 

-0.755 

- 8.055 

0 

0 


v, 

*£ 


0 

-19.4 

-1376 

-12A 

-299.4 

- 0. 755 

-8.055 

0 



oC 


55.11 

-1375 

-60.62 

-299.8 

-14.33 

- 8.005 

-.177 

0 


*3 

« ; 


83.32 

-Z057 

29.74 

8946 

16.95 

283.4 

0 

0 



i 


0 

83.32 

-2057 

29.74 

8946 

16.95 

283.4- 

0 


% 

\ 


10.63 

-262-4 

19.42 

2140 

1 2.9 

407.6 

0 

0 





0 

10.6 3 

-261.5 

19.42 

2140 

12 . 9 

407.6 

0 



& 

/ 


9.014 

-224 

33.23 

3008 

16.66 

664.1 

1. &16> 

17.9 


*8 


(4-21) 


Model 

The model to be used with this example will be defined by the third- 
order system 


*» 


0 

1 0 




= 

0 

0 1 


n 

A. 


0 

-1 -14 \ 


^3 


(4-22) 


The model incorporates the concept of the drift minimum require- 
ment, the condition that leads to the cancellation of the sum of all force com- 
ponents, such as gust inputs, perpendicular to the nominal flight plane. 

This cancellation of forces is equivalent to requiring a free integrator in the 
closed-loop characteristic polynomial. Therefore the model characteristic 
equation contains a root at the origin of the s plane 

| Js- Lj = s(s 2 + 1.4 s + t). 

Analysis of the Optimal System 


Substituting the characteristic polynomials of the open-loop booster 
and the model along with their adjoints in Equation 4-13, the root square 
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locus expression, the result is 

0=i .l s 2 (s 2 * 1.4 s* f) (4-23) 

r ( Si 17. $)(s t . 2 94)(s i . in) (s t. 04ns) [s 1 tz(. 005)z. 32 54 (2, 32)*} [s 2 ^ Z(,00S)5M S 4 (.5.64) Z ] 

The locus of roots of the optimal system is plotted in Figure 4.1. The 
adjoint, or right-half plane part of the locus is omitted for clarity. The lo- 
cus shows that the open-loop roots of the booster originally associated with 
the rigid-body poles tend toward the model singularities, shown in the locus 
as zeros. The remaining poles, originally associated with the bending modes, 
tend to become distributed in a Butterworth filter pattern as the parameter 
of the locus, q/r, becomes large. For any value of the weighting of the 
error portion of the integral to the control, the accuracy of the approxima- 
tion to the model in terms of the closed-loop root locations can be determined. 


Optimal Control Law 


The optimal control law is of the form /3 C = - !<?% where K is a 1 x 8 
matrix of feedback gains from the state variables of Equation 4-20 to the 
control inputs. The closed-loop characteristic polynomial is 

A(s) = lls-F + G>c\ 

~ 5 n * (b„ r £ n ) s* r + (b n .f + ^rj-f) 5* ' (b>2 6 'fc-z) S + ( b, + £. t ) 

» S 8 + (I7.99*t e )s 7 + {38.71 4 (M6'.S + £js s + (f93.7 + & s )s 4 (4-24) 

+ (3010* 4-(20.7t-h£. s )s Z 4 (-225.5' + djs + (9.074 + 


The closed-loop characteristic polynomial can also be obtained direct- 
ly from the root square locus plot of Figure 4.1 for any value of q/r. For 
instance, assuming that q/ r = 100 will yield an optimal system that results 
in an acceptable approximation to the model, the closed-loop characteristic 
polynomial is 

A(s) = (s+.0077)(s4/7.9)[s 2 -i- 2(.7s)(.82)S4-(.82) 2 '] [s 2 4 Z(46)(2.6)s *(2.d) Z '\ [s 2 / 2(./2)(^2) S + ^.z) 2 ] 

, (4-25) 

= S 8 f 23.02 s 7 + 145.87s 6 + 1135.8 s 5 4 3372.5s 4 4 788S.6s 3 + 7157-4 s z + 3107.8 s * 24.11 

Equating coefficients of the powers of s in Equations 4-24 and 4-25 yields 
the feedback control law 

/3 t = -15.0m, - 3331.3 % z - 7138. 7%y 4875.6% - 3178.7 V s ~ 470.3M u - 1072 ^- 5 . 03 ^ (4-26) 

Or, in terms of the variables of Equation 4-21, the control law 
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igure 4. 1 Locus of Poles of Optimal System (Drift Minimum Model) 
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becomes 

$=-£V = -^/rV 

« 1.93Z0 e f S.3G0 e ^.893oc- -239^--^2 i li ~.^SGri 2 -.042lj i - .Z&1 (3 (4-27) 

With the exception of {3 the variables of Equation 4-27 are still not in 
terms of directly measurable quantities, i. e. , in quantities that directly 
represent outputs of sensors located on the body of the vehicle. A set of 
measurable quantities can be chosen, such as the vector 

^ 5 ^ 7 ^ 1 ^ ’ ^-1 (4 - 28 ) 

where 0 %1 , 0 y, , f and , represent outputs of a position gyro, 

a rate gyro and a normal accelerometer located at station % f on the 
vehicle body and 

A0 1 - 0^ - ^ 

A 0 2 = 0 Vi - 0^ 

A 0, = 0 %i - 0 %z 


A transformation = 13y can be formulated that defines the meas- 
urable quantities of Equation 4-28 in terms of the variables of Equation 
4-27. With this transformation, the feedback control law becomes 


/3 0 = - k?A~’ 

If one chooses the vehicle body stations 


(4-29) 


Vi - 41. 5 meters 

v i = 86 meters 




122.5 meters 


the feedback control law becomes 

fa = 1.933 +&S6 0P + 0.4-24 a ^ + 72.98^0^0.208 I9MA0> Z + -0329 A0> z (4-30) 

' ' ** +3.06(3 

The differential gyro realization of the feedback control law for this 
particular case results in high feedback gains from several of the measured 
quantities, resulting in a feedback control law that may be difficult to mech- 
anize. It would be more realistic to select another set of measurable quan- 
tities such as the set 


% - [(**, • • v, H< K • ) ( K . v, ) A 


(4-31) 
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If this set is chosen, the control law becomes 
(3 C = f.93J 0^ +£48 ^4.0466 -.6&4 0 %g +M96&^ 


sensors located at sensors located at 
stat i on y f - iH . 5m stat ion % x = 86. Om 


where the units of (3 / 0 and / 0 are rad/rad and rad/ (rad/sec) and the 
units of /3 / CLy are rad/ (meter /sec 2 ). The control law of Equation 4-32 in- 
volves feedback gains more easily mechanized than those of Equation 4-30, 
but either control law will, of course, result in the same closed-loop char- 
acteristic polynomial and closed-loop transient response. It is clear that 
the realization of the optimal system is not unique, and many engineering 
factors enter into the choice of sensed quantities and sensor locations. 


- * 0592 0^ + .0244 a ^ - 0. 5S/3 l 




sensors located at 
station y = 122.5m 


(4-32) 


Equivalent Compensation Network 


The optimal feedback control law may be synthesized as a filter or 
compensation network in the feedback path which shapes the signal of the 
output of one or a linear combination of sensors. There are, however, def- 
inite restrictions and qualifications that must be placed upon such a synthesis 
procedure. The two systems, one synthesized by feedback gains and the other 
synthesized through a filter network, will not respond exactly alike to initial 
conditions of the state. Therefore, one must be content, when using a feed- 
back filter synthesis procedure, with a synthesis of the closed-loop poles of 
the optimal system rather than a complete time history equivalence. 

A filter must satisfy certain characteristics of realizability that simply 
do not have to be considered when gains are fed back from a set of state var- 
iables. So the filter problem is more difficult than the straightforward feed- 
back gain procedure of synthesis. The technique of specifying a filter is 
uncomplicated. Because the system has been assumed to be completely lin- 
ear, the state variables are related to each other through transfer functions 
and therefore one state variable can be reconstructed from a measurement 
of a second state through a filter network (with the exception, of course, of 
the initial transients of the filter itself)- There are limitations to this kind 
of reconstruction. The most important limitation is that the open-loop trans- 
fer function of the measured state variable must be minimum phase. Other- 
wise, the filter will have right-half plane poles, creating an unstable config- 
uration for even the slightest variation of the launch vehicle or compensation 
network parameters. The second limitation to the reconstruction procedure 
is that the transfer function of the measured state should have a numerator 
polynomial of order n - 1, where n is the number of first-order equations of 
motion that describe the system's dynamic motions. This requirement will 
quarantee that the compensation network will contain a numerator polynomial 
that is not of higher order than the denominator polynomial and the filter will 
be realizable. 
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hi 

Each of the transfer functions (s) of the open-loop launch vehicle 
can be expressed as the ratio of two polynomials in s of the form 

fie (,) ' V <S> 

where 


( s ) = the numerator polynomial of the transfer function for 
the ith state variable, and 

^ ( 5 ) = the open-loop characteristic polynomial of the launch 

vehicle. 

The transfer function of the state variable is simply related to the trans- 
fer function of the state variable through the expression 





(4-33) 


Using this relationship among the states, a feedback control law 
A = can be reconstructed from a measurement of a single state, say 

, in the following manner 








^ + ^ N (s) 3, * ^3 ~r ^ % + * • • + —y- ( s ) 'ff 


(4-34) 




*1 




= ~ (S) X ^ ^ 

'Yf * = / 


The transfer function of the required compensation network is then 
given by 


Y(*) = 


t 1 f 


1 

A/, (s) 


L A (*) 

i=Z 


(4-35) 


Suppose it is desirable to reconstruct the control law in the manner 
described above from a measurement of not more than three of the state 
variables of Equation 4-32. Consider the numerators of the transfer 
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functions of these states which are given below: 


'Xi 


(*,) 


A U 




Nk 


v* 2 ) 


/X 




Na . 


(*■*) 


N a 




(s) = -19.4 s 5 - /. 335 s 4 - 5.70 s 3 - 13-6s 2 - 1315s - 21-85 

(s) - -19.4 s 6 - 1.335s 5 - 5.70s 4 - 136 s 3 - 1515 s 2 - 21-855 

(s) = - 141.1s 7 + 1794s* 1 2131 S 5 4 164. 3 S* 4 28,260 s 3 - 33.70 s 2 - 9450 s * 382.8 

(s) - Z 3. 49s 6 4. 44 3s 5 4 52.65s 4 - 10.69s 3 - 1443s 2 - 16.68s 

(s) = 62.65s 7 - 16.46s 6 - 11,460s S - 4H.95 4 - 45,475 s 3 - 966- 9s 2 - 6 774s 4 274 .6 

(s) = - 31.48 s 6 -. 5/46 s 5 4 494.9 s* -11.72 s 3 -1591s 2 -10. 67 s 

(s) = 255.7s 7 - 6 . 723 s 6 4 1764s 5 - 784.3s *- 110,406s 3 -1431 S 2 - 364/ s + 147.7 


(4-36) 


V* (s) = i7.9s 7 4 1.616s 6 4 664.1 s S iZU6s* 43008s 3 4 33.23s 2 - 224s 4 9.074 


We wish to choose as our measured quantity a state or combination 
of states whose transfer function contains a minimum phase numerator 
polynomial of not more than seventh order. None of the polynomials of Equa- 
tion 4-36 qualify in themselves, but a combination of these polynomials will 
satisfy the measurability requirements. Several combinations will yield the 
desired characteristics, so as an example we can choose the sensed quantity 
to be A = fi - ( <2 < > Xl 4 <2> Vi ) . The numerator of the transfer function 
is given by 

(s) = Np (s) - ( s) t Nfy' j fo)] (4-37) 

Substituting for (s), {s) and A/<ji^ (s) from Equations 4-36 

A/ 4 (s) = 17.9 s 7 + 21.02 s 6 + 68 S. s *4 598 s* 4 3590 s *4 / 362 s*+ 1 / 13 S + 30.9 ( 4 - 38 ) 

N a (s) is of seventh order and a Routh's criteria check of N » (s) =■ 0 shows 
all the roots to have negative real parts. 

Using this sensed quantity, there are still several ways the feedback 
control law may be synthesized. For instance the quantity A - + 

may be the only quantity fed back to the gimbal actuator input, as shown in 
Figure 4. 2. 
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Figure 4. 2 One Realization of the Control Law 

or, as a second example, the attitude and attitude rate parts of the control 
law may be independently fed back, as in Figure 4. 3 below: 



Figure 4.3 A Second Realization of the Control Law 


Using the results of Equation 4-34, the transfer function of the com- 
pensation network of Figure 4. 2, which realizes the control law of Equation 
4-32 is given by Equation 4-39: 
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Witt) 1 


T 

1--1 M a (s) N A (s) 


/• 933 N 6 (s ) * U8 Nq ( 5 ) * .0486 ty, (s)-,0S6+ty (s) - . 01% N a ($) 

{% 1* (%1> *(%i) ( *l) Wj) 


-.O59ZH0 (s) + • 0244- N 0 (s) - 0. 35 N* (s) 

( *s> f(y s ) 


-5.30s 7 - 106,8 s 6 - 469-Zs*- 3168.8 s 4 - 4880. 3 s 3 - 7/gg.Ss 2 - 33Z7s - 15.02 
17.9s 7 * 21.02 S 6 * 685. S S * 598s 4 * 3590 s 3 * /362s z *///3S *■ 30.9 


(4-39) 


-530(s*. 0456, )(s*. 635)( 5+ 16. 8 1) [s Z *2(. 163)(5. 196) S* (5.)%) Z ][$ 2 * Z( 33?)(l.468 ) 5* (f.468) 1 ] 


f7.9(s*.0287)js Z * 2(.0219)(5.64-)s*(5.64) z ][s 2 * 2(fl4-)(2M)s*( 2. 3d) 2 ][s Z * Z(.308)(.577)s + (S77) z ] 


The transfer function of the compensation network for the realization 
shown in Figure 4. 3, on the other hand, is given by 


8 


g _ /_ 

i-5 % (S) % (S) 


.0486 h!„ ( 5 ) --05641/A. (s)*.0Z%Ha (s)-.0S92hlk ( 5 ) *-0244/1, (s)-ASSAi(s) 

f(% ) ‘ X /) If* A (x.\ fa.. \ ^ 










- 6.30 s 7 - ,587s i ’- 424.4s S - 46.2 s f - 3705. 7s 3 * 97.84 S Z - 667.4- 5 + 27,18 
17.9s 7 * 21.02 s 7, *685s s *598s*+ 3590s 3 * 1362s 2 * 1f/3s *-30.9 


(4-40) 


-5.30(s-.04/5)[s Z *2(-.0O12j(8.37)s*(a37) Z ][s Z *2f.0256)(3./25)s*f3.12s) Z ][ 5 \2f0/45)(4296)s*f42%) 2 ] 
17.9 (s* ,0Z87)[s* Z(. 0W)(5M) 5* (5.64) Z ][s Z * Z(. m)(l 38)s + (Z 38) 2 ][s 2 * Z(. 308) (.577) 5 * (■ 577) Z 1 


Although the filter transfer functions described in Equations 4-39 and 
4-40 are of seventh order, they may be realized without great difficulty. If 
it is objectionable to have a seventh-order numerator polynomial, then it 
may be possible to synthesize Equation 4-30 instead of Equation 4-32 using 
the same sensed quantity A = /& - * 0*/) but feeding back the actuator 

deflection part of the control law 3. 06/3 separately (i. e. , not reconstructing 
its effects by the compensation network). The resulting compensation net- 
work transfer function would be the ratio of a sixth- to a seventh-order 
polynomial. 
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So again it is found that the synthesis of the control law is not unique, 
and a search should be made of the control configuration that presents the 
fewest problems from a mechanization point of view. 

The compensation network transfer functions of Equations 4-29 and 
4-30 were synthesized from expressions of the exact optimal control law of 
Equation 4-32. If a compensation were to be synthesized for an actual 
vehicle, there are several steps that would be taken in the normal event of 
system design that are not shown in the preceding analysis. First, approx- 
imations to the control law of Equation 4-32 would be made. Perhaps some 
of the feedback gains can be eliminated entirely, a procedure that would alter 
the strict optimal character of the design, but may still yield a satisfactory 
and acceptable vehicle response. The second procedure would be to simplify 
the compensation network itself to yield a simpler design, still retaining 
an acceptable closed-loop vehicle response. It is felt that a fair number of 
simplifications can be made before the response becomes unacceptable (such 
as unstable bending modes). Therefore, the optimal control law can be used 
as a guide which defines a superior system. The aim is to realize a control 
law as close to the optimum as practical considerations will allow. 


Multiple-F eedback, Multiple-Compensation Synthesi s 


Although it may be possible to synthesize the control either as feed- 
back gains or as a single compensation network, neither of the two synthesis 
techniques are entirely satisfactory. The first technique requires a feedback 
path for each state variable used to describe the dynamic characteristics of 
the launch vehicle. It is clear that when many elastic modes of motion are 
to be controlled, the number of sensors required to realize the synthesis is 
prohibitive. In addition, there would be no way to compensate for the dy- 
namics of the sensors themselves, which may present a problem if one is 
trying to control bending modes whose natural frequencies are in the neighbor- 
hood of the natural frequencies of the sensors themselves. Approximations 
to the control law, obtained by ignoring feedback paths or by using least 
squares approximations, may still require too many feedback paths or may 
result in closed-loop dynamics that are altered far beyond allowable toler- 
ances. 


The second technique, which requires that a compensation network be 
synthesized to shape the signal output of a single sensor, is an equally unsatis- 
factory approach to optimal control system synthesis. If the compensation 
network is synthesized in the feedforward path with unity feedback, the com- 
pensation network must contain zeros that cancel the unstable open-loop poles 
of the launch vehicle. Because pole-zero cancellation of unstable poles is phy- 
sically impossible, the technique is worth not more than ephemeral consider- 
ation. No systematic approximations to the optimal will produce a mechaniz- 
able, stable result without at least an increase in the order of the system. 

A compensation network in the feedback path shaping the output of only one 
sensor has difficulties almost as insurmountable as the feedforward compensa- 
tion scheme. If the transfer function associated with the measured output is 
non-minimum phase, the exact compensation network must contain right-half 
plane poles to cancel the right-half plane zeros of the transfer function of the 


42 



IH-2089-F-1 


sensed quantity. Therefore, in order to realize this type of compensation, 
one must either choose a sensed quantity whose transfer function is mini- 
mum phase or one must approximate the filter, either by deleting the unstable 
portion of the network or by approximating the unstable part with stable com- 
ponents, a procedure that may produce some difficulties. In any case, the 
order of the filter before approximations are made is equal to the number of 
state variables used to describe the launch vehicle dynamics. If a dozen 
elastic and slosh modes of motion are to be stabilized, the order of the re- 
quired filter can be very high. Finally, the stability of the launch vehicle 
depends upon the reliability of sensing a single dynamic quantity, which may 
not be desirable. 

A review of the difficulties associated with the synthesis procedures 
outlined above can lead one to the requirements for a desirable synthesis 
procedure: 

1. The synthesis procedure should involve more than one sensor 
but fewer sensors than the number of states of the system. 

2. The compensation networks should be in the feedback paths 
to avoid the pole-zero cancellation problems associated with 
feedforward compensation, 

3. Each compensation network should be of order lower than 
n-1, where n - the order of the plant. 

4. The synthesis procedure should be methodical, with systematic 
selection of the poles and zeros of the compensation networks. 

5. The synthesis procedure should contain the ability to compen- 
sate for sensor dynamics. 

6. The synthesis procedure should be amenable to rational, 
systematic methods of simplification. It should be possible 
to predict beforehand the consequences of a particular step 
in the simplification process. 

7. The order of the characteristic polynomial of the closed-loop 
system should theoretically equal the order of the open-loop 
system, even after approximations and simplifications are made 
to the compensation networks. 

8. The compensation networks should be chosen, if possible, such 
that the vehicle would not be unstable with the failure of any one 
feedback path. 

The requirements dictate that the synthesis procedure should be a 
compromise between the feedback gain procedure, which avoids the pole- 
zero cancellation problem, and the single sensor synthesis technique, which 
eliminates the need for an excess number of feedback paths. Schematically, 
the block diagram of Figure 4.4 indicates the requirements: 
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Figure 4. 4 Block Diagram of Desirable Control Law 


where the sum of the poles of the three compensation networks equals n - 3, 
where n is the order of the system. 

The requirement that a systematic method be devised for the selec- 
tion of the filter poles and the desire not to increase the order of the closed- 
loop system demands that the poles of the compensation networks be chosen 
from among the zeros of the C s ) transfer functions. 

The dynamical description of the elastic launch vehicle must be such 
that three of the states are , ^3 ■ The remaining states must be 

chosen such that each state is related to y t , , or by compensation net- 

work component parts. In other words, we wish to be able to describe the 
control law as 


a = *;yi 

* " - ••• V* **] 
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The control law can then be specified in the form; 
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The control law of Equation 4-43 specifies three compensation networks 
shaping the outputs of the three sensed quantities y, , y t , and y 3 . Because the 
order of the filters is p < n - 1, q < n - 1, r < n - 1; p + q + r = n-3, the 
poles of the filters may be selectively chosen from the roots of ( S * ), 

( 5 + /3j ), and ( s + ^ ), which have negative real parts, avoiding the pole -zero 
cancellation problems of single - sensor synthesis. 


The control law of Equation 4-43 suggests a methodical sequence one 
may use in the process of simplifying or approximating the optimal control 
law. The contribution of each state variable to the total compensation net- 
work depends upon the magnitude of the feedback gains associated with the 
particular state. For instance, if in the compensation network 


Y (s) 


-i z &, c 

+ + 

S+'flfcy S + 


+ 'jP^± 

S + CCp 


(4-44) 


the feedback gain k ? is found to add negligibly to the closed-loop dynamics, 
it may be possible to ignore the gain, yielding the approximate network 


resulting in a simpler, 


f S + 06 2 S + QCp 

lower order compensation network. 


(4-45 ) 


Required Dynamical Description of the Launch Vehicle 


In order to synthesize the system with three measured quantities and 
three compensation networks, it is probably necessary, and certainly exped- 
ient, to describe the system in the following form: 
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(4-47) 


The elements of the F matrix are obtained from the open-loop char- 
acteristic polynomial D(s) = + . *■ d 0 . The elements of 

the A matrix are obtained as follows: The first row of A is obtained from 

the numerator of the (s) transfer function Q, v , f s n ' 1 *(Z n . l s”' 2 + • • - • 

The second row of A is obtained from 
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Similarly, the third row is obtained from 
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and all other rows of A are obtained in a similar manner. 


Simple Example 


Consider the plant represented by the block diagram 



It is desired to synthesize the system as shown, with two feedback 
paths and with a realizable compensation network in the y 2 feedback path, 
such that the closed-loop characteristic polynomial is 

A( s) = s 3 + 2s 1 + 3 s + 4. 


From the transfer functions 



S z -a s - 2 
S 3 -h s l + 5 + 2 


y z s z *■ 7s *-12 

— (s) = — 

Li S 3 + S z + 5 + 2 


the state space description of the plant is obtained 
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where the first row of A is obtained from numerator of the fy/a (s j transfer 
function s 2 + s - 2 and the second row is obtained from 


S z + 5 - Z 
s + z 


= s - 1 


while the third row is obtained from the numerator of the (s) transfer 
function s 2 + 7s + 12. 

Solving for u = -Kx, we follow the procedure of Equation 4-24 


-2 - kj = -4 

-1 - 4 = -3 

-1 - k 3 - -2 


V 2 

k 2= 2 
k 3= 1 


In terms of the variables ^ , y 1 f , and resu ^ ls 
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Since 


V, ( (*) - t/ 9+ 2 , 

u-(s) = - 


the control may be synthesized as 


i]z_ 

S + Z 


<j, & - V ^ 


% 5 ^/ 

L 5^-2 


0 , (s) +■ — y z (s) 

4 


which yields the desired closed-loop characteristic polynomial 

A (s ) = s ^ + Zs 2 + 3s + 4. 


Although the example is quite simple, it does illustrate the princi- 
ples involved. The number of feedback paths is less than the number of 
states, yet greater than one. A compensation network has been systemati- 
cally specified for the non-minimum phase variable y^, yet the compensation 
network is realizable, containing no right-half plane poles. In addition, it 
would be relatively easy to investigate the consequences of omitting a feed- 
back gain, but a simplification of this type would not increase the order of 
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the closed-loop system. 

The multiple -feedback, multiple -network synthesis procedure was not 
tried in the elastic launch vehicle problem, but it is clear that the techniques 
illustrated are directly applicable. 


Response of the Optimal System 


The closed-loop response of the vehicle, using the feedback control 
laws of Equation 4-32 or Equation 4-30, is shown in Figures 4.5, 4.6 and 
4.7. Figure 4.5 shows the normal acceleration response at the three sensor 
locations indicated in Equation 4-32. The input command was a step, a 
severe excitation to such a system, but calculated to excite the bending mo- 
tion and enhance the residues of the poles originally associated with the 
bending modes. The effect of the first bending mode shows clearly only at 
station % 5 = 122.5 meters. The second mode oscillations are clearly shown 
at all three locations. These oscillations are reasonably well damped (for 
a bending mode), the measured damping ratio of this mode closely matches 
the predicted X = • 1 2 from the root square locus plot. Figure 4.6 shows 
the attitude response of the vehicle at the c. g. , % - 41.5 meters. The re- 
sponse is smooth, showing very little contribution from either bending mode, 
but this is to be expected because the mode slopes are small at this location. 
Figure 4.7 shows the control motions required to achieve the responses 
shown in the previous two figures. Except for the initial abrupt response, 
the control motion is reasonable, and shows, in general, that the system 
could probably have been speeded up somewhat before encountering an un- 
acceptable control motion. 


Parameter Variations 


The bending characteristics of a large flexible launch vehicle are 
known only to a fair degree of accuracy. The normal bending mode shapes 
and slopes are subject to uncertainties depending upon the mode number and 
the vehicle body station. It is necessary, therefore, to calculate the effects 
of variations of the elastic properties of the vehicle on the closed-loop opti- 
mal system. A series of variations in these elastic properties was speci- 
fied and their effect on the system, defined by the control law of Equation 
4-30, was calculated. The transformation B relating the vector y to the 
vector tj is defined by the equation shown at the top of page 53. 
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Figure 4.6 Attitude Response of Optimal System 
at Booster Center of Gravity 



Figure 4. 7 Time History of Optimal Control Motion 
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where 


Y: 




t h th 

normalized mode deflection of the j mode at the i body 
station 


^ th th 

Vj ^ = normalized mode slope of the j mode at the i body station 

At body stations x 1 = 41. 5 in. , X z = 86. 0 in. and k 3 = 122, 5 in. , the 
nominal bending mode and slope values and their selected variations from 
the nominal are given in Table 4. 1 below: 


TABLE 4. 1 

NOMINAL MODE CHARACTERISTICS AND 
EXPECTED RANGE OF VARIATIONS 

Mode Characteristic Nominal Value Expected Range of Variation 



-0. 390 

±15% 


+ 0. 028 

±15% 


-0. 555 

±15% 

y' 

Y ^V 

+0. 0085 

±15% 

Y' 

Y< ^ 4 ) 

-0. 025 

±25% 

y' 

-0. 060 

±35% 

y' 

T U% t ) 

-0. 123 

±25% 

y' 

+ 0. 143 

±25% 
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A brief program was conducted to investigate the variations of the 
closed-loop poles of the system with the control law fixed as defined by Equa- 
tion 4-30. The results for a few of the parameter variations are tabulated in 
Table 4. 2 below: 


TABLE 4. 2 

VARIATIONS OF THE CLOSED-LOOP POLES 


NOMINAL 

ALL MODE 
PARAMETERS 
INCREASED 
(SEE TABLE 
4.1) 

ALL MODE 
PARAMETERS 
DECREASED 
(SEE TABLE 
4. 1) 

Y f 

Y *(x,) 

ONLY 

INCREASED 

V' Y' 

INCREASED 

s b -.00770 

-.00833 

-.00790 

-.00750 

-.00820 

s = -17.901 

-17.901 

-17.901 

-17.901 

-17.901 

s = -.612 ± j .547 

-.614 ± j.478 

-.615 ± j.526 

-.611 ± j - 565 

-.618 ± j .496 

s = -1.200 ± j2.304 

-1. 187 ± j 2. 319 

-1. 197 ± j2. 309 

-1.200 ± j 2. 299 

-1. 193 ± J 2. 3 16 

s = -.743 ± j6. 155 

-.744 ± j6. 154 

-.744 ± j6. 155 

-.744 ± j6. 156 

-.744 ± j6. 1515 

A casual examination of Table 4. 2 will disclose that the optimal system 
is not particularly sensitive to variations of the bending mode parameters of 
the launch vehicle, even though the variations of the closed-loop roots tabulated 
in Table 4.2 represent the widest change of closed-loop dynamics discovered 
using the selected mode characteristics of Table 4. 1. The variations were 


selected at random, and therefore an absolute conclusion concerning insensi- 
tivity cannot be made. However, all evidence points to a very insensitive 
de sign. 


Control Law Approximation 


A least squares fit of the control law of Equation 4-30 using only three 
of the eight state variables was calculated to determine if a direct approxi- 
mation to the exact control law could result in an acceptable elastic booster 
control system design. 

The requirement is that the outputs of three sensors located at station 
=41.5 meters are to be used to approximate the vector 

</' = [A . ^ > h • 4] 

and therefore the optimal control law of Equation 4-30. The relationships 
defining the outputs of these three sensors are: 
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(4-49) 


where B. consists of the first three rows of the square matrix B defined by 
Equation 4-33. 

The least squares fit is calculated by computing the generalized in- 
verse of B^ called B^ (Ref. 7). The generalized inverse is obtained by 
solving for the matrices N and M such that 

, , , (4-50) 

S/s; M =N73,B,' = 3/ 

then 

3* = NB 1 M 

The approximated feedback control law then becomes 

8 . (4-si) 

For this particular problem, the least squares fit to the control law 
of Equation 4-30 is 


f3 a = 1. 93 95 0^ + 6. 302 0 %i ~ 0. 0103$ (4-52) 

The feedback control law of Equation 4-52 can be easily mechanized. 
The closed-loop characteristic polynomial obtained by using the approximate 
optimal control law of Equation 4-33 is; 

A(s) = ( 5 + . 02 Z$j(s+. 377 )[_S 4 .428 ±j/.56][s + 6.95 + j4.6d] [s + .703± j 5.oj] 

Comparing this characteristic polynomial with the characteristic 
polynomial of Equation 4-21 shows that the approximation has preserved sta- 
bility but the drift minimum property of the design has been compromised. 

It appears, however, that the technique shows promise for obtaining accept- 
able control system design approximations of the optimal system, although 
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there is no reason to believe that all approximations of the type given above 
will result in stable closed-loop configurations. 

Perhaps the best way to incorporate the idea of the generalized in- 
verse into the control problem is to use the generalized inverse theory to 
adjust the zeros of the compensation networks of Equation 4-45 after the 
k. having negligible effect on the closed loop dynamics have been deleted from 
tlie control law. 


Addition of the Third Mode 


A third bending mode was added to the definition of the launch vehicle 
equations of motion. The purpose of this addition was to obtain the effect on 
the third mode of a control system designed to increase the damping ratios 
of only the first two modes. The equations, which include the third mode, 
are given in Appendix A. The third mode has a natural frequency of & = 9. 18 
rad/ sec and a damping ratio £ = . 005. Using the control law of Equation 4-28, 
the closed-loop poles originating from the third mode were unstable, having 
a closed-loop natural frequency of = 7.885 and a damping ratio of X> - -.206. 
This mode could, of course, be included in the system design and stabilized 
with a corresponding increase in the complexity of the control law. 

It should be stated that the sensor positions selected for the control 
law of Equation 4-32 were chosen without regard to the third bending mode 
characteristics at the particular sensor locations. An obvious way of re- 
ducing the effect of the control law on higher modes is to judiciously locate 
the sensors to minimize the mode pickup. A less obvious way is to add a 
constraint to the problem that would have as an objective a design that would 
decrease the sensitivity of the third mode. This design would result in a 
filter to desensitize the third mode. 

A few generalities can be made on the use of a model to specify a 
desirable response. The technique seems to be the straightforward, sim- 
plified way to design a closed-loop system, optimal or otherwise. The 
optimal aspect of the problem is secondary if the only requirement is to 
match the poles of the closed-loop system to the poles of a model. The opti- 
mal aspect of the selection of approximations to the model should be pre- 
served for at least three reasons which are implied, but not proven, in this 
section. First, the optimal formulation of the problem considers control 
motion deflections. If the control motion amplitudes are too large (for a 
specific input) it is necessary only to weight the control portion of the per- 
formance index more heavily, resulting in lower amplitude control motions. 

No other design technique, particularly for complex systems, can specifi- 
cally consider the control motions as part of the problem objectives. A 
second advantage of optimal control is also peculiar to the technique. It 
has been observed that the closed-loop poles rapidly approach the Butter- 
worth configuration, i. e. , the asymptotes of the locus. In general, the 
root square locus plot appears to approach its asymptotic values faster than 
for an ordinary root locus plot. Once the closed-loop poles approach the 
asymptotic values, the poles become relatively insensitive to changes of the 
coefficients of the original equations of motion. Since a large number of 
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examples have shown that the root square locus approaches its asymptotes 
at a lower increment in frequency than an equivalent root locus, it can be 
safely stated that a linear optimal system is equivalently less sensitive than 
a conventionally designed system. Finally, the root square locus expression 
enables the designer to systematically select all of the poles of the system, 
particularly when the model and the plant are of different order. 
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SECTION 5 

THE APPLICATION OF LINEAR OPTIMAL CONTROL 
TO SYSTEMS CONTAINING UNCERTAIN PARAMETERS 


Introduction 


Most of the synthesis techniques available to the designer depend 
upon the assumption that the plant parameters are well known. This assump- 
tion is not valid when one considers the design of a controller for a flexible 
booster. Variations in the plant parameters are, of course, expected to 
occur over the flight envelope of the booster, but quite aside from this effect, 
there is the problem that even at a specific flight condition (where one can 
usually assume the plant parameters are constant) the values of the various 
coefficients, particularly the flexible mode shapes and slopes, are not known 
exactly. 


The bulk of the theory which treats this case has, in the past, consi- 
dered the situation in which the statistical variations of the plant parameters 
take place over a period of time. The situation which is of more interest is 
the case where the plant parameters are relatively constant over time but 
can be considered as random variables described on a sample space. Some 
recent research in this area has approached the problem by incorporating 
some measure of the randomness or variability into the performance index 
prior to the minimization and specification of an optimal control law. In this 
way an expression for the optimal control is developed which reflects the 
fact that some of the plant parameters may not be at their nominal value. 

One approach to this problem can be formulated by asking for the 
compensating networks which minimize the E { 2 Vj. That is, 

r 

J ( tj f Q L/ + u, ' 13 u )dt 

0 

If the equations of motion which describe the system are manipulated 
into the vector form 

% - F % + Qti' 


win. f . min- 

conip* E [ZVj = comp- E 

n etu>*rk netujork 



where F and G have entries which are random variables with the joint density 
p(F, G), this problem can be stated as 


min. 

comp. 

netutork 


r 



dfdG 


In this section, it is not our intent to formulate the problem in the most 
general manner possible (that is, consider the multi-controller and multi- 
output situation) since it is desirable that we gradually work our way into the 
problem and develop a "feel 1 for the mathematics involved. It will suffice to 


*In this report expressions such as E ( 2 V J will read as "The control 
u which minimizes the expected value of 2y. 
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say that all the tools required to extend the single-input, single-output case 
to the multi-controller case are already available within the present frame- 
work of the calculus of variations. Thus we will be content with outlining 
the background philosophy, solving three simple illustrative problems, and 
applying what has been developed to the task of investigating the character- 
istics of the optimum compensating network required for a flexible booster 
(only the first bending mode will be considered). 

The equations which define the conditions of optimality will be de- 
rived using frequency domain methods. Appendix B considers some of the 
difficult points which arise with nonminimum phase systems since the trans- 
fer functions which describe the booster have both poles and zeros in the 
right-half plane. 


Preliminary Comments 


To begin, assume the feedback control system shown in Figure 5. 1. 


16 


<? 


VJr 


U 


W 


Figure 5. 1 Forward Loop Compensation 

The compensating network, W , is to be designed in such a fashion that the 
expression 

oo 

0 ' 

is minimized. In Figure 5.1, Lc represents the control (or input) to the 
fixed elements of the system W(s). The above statement may be expressed 
in the form 

0 

Although it may be possible, it is believed that the analytical difficulties 
involved in minimizing Equation 5-1 by solving directly for the optimum 
compensating network are insurmountable. In the past, some writers (Refer- 
ence 8) have attempted to sidestep the difficulties involved by reformulating 
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the problem in an open-loop manner, and then solving for an open-loop com- 
pensating network. That is, they formulated the problem as depicted in 
Figure 5.2. 



Figure 5.2 Open- Loop Configuration 


This, of course, is completely unsatisfactory when W(s) contains nonminimum 
phase components. 

In the event that W(s) is completely known, the analytical difficulties 
involved in finding the optimum compensating network may be avoided by 
solving for the optimum value of some other system variable such as the 
optimum control or the optimal error. These problems turn out to be quite 
tractable and present no insurmountable analytical road-blocks. Once either 
the optimal error or the optimal control has been found, it is a simple matter 
of algebra to solveTor the optimal compensating network. As an alternate 
procedure, one might elect to solve for both the optimal error and optimal 
control and then define 


oo 

f l (e‘*ru.’)dt 

f ( e2> ru. z )M 

0 

All of these procedures obviously yield the same expression for W when W 
is completely deterministic. 

When W(s) contains uncertain coefficients, described only by proba- 
bility density functions, it seems reasonable to find the compensating network 
which minimizes the expected value of the performance index. However, 
one cannot solve for just the optimal control which minimizes the E {2V} and 
then algebraically manipulate the blocks to find W since W is unknown (for 
any given system). Nevertheless, it can still ma£e sense to define 


K = 


s[j g (e 2 +ru z )dt^ 

(e z +■ ru> z ) dt 


(5-3) 


That is, define the compensating network to be the ratio of the control which 
minimizes the expected value of the performance index to the error which 
minimizes the expected value of the performance index. Defining the 
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compensation in this manner is a subjective decision which is certainly open 
to discussion. However, one feels that the final decision as to the "goodness" 
of the choice should rest with the results achieved — namely, does it design 
"acceptable” systems. 

It is, of course, difficult enough to decide when a system is accept- 
able even when all the parameters are known precisely — thus we expect 
that our difficulty in deciding this question must increase as our knowledge 
of the system becomes less precise. However, the simple illustrative exam- 
ples which follow do seem to indicate that there is merit in the idea of finding 
the compensating network which minimizes the expected value of the perfor- 
mance index. Notice that this does not necessarily imply that the perfor- 
mance will be less sensitive to parameter variations nor does it guarantee 
that any particular system which uses the average compensation prescribed 
by the analysis will be stable. (Of course, one would hope that this might be 
a by-product of the design procedure. ) It would thus seem worthwhile to 
expend additional effort to determine the relationships, if any, which exist 
between parameter sensitivity, stability and the minimum of the expected 
value of the performance index. 

We will take Equation 5-3 as the definition of the compensating net- 
work to be employed for the system depicted in Figure 5. 1. If a different 
feedback configuration is desired, say that of Figure 5. 3, W might be defined 
as c 

12- E{f“(e z +ru*)dt} 

K- T E {<C(e*<-ru,*)dt} (5 _ 4) 

where 



Figure 5. 3 Feedback Compensation 


It is still necessary, in any alternate compensating scheme, to solve 
for at least two optimal variables of the system. 

In the deterministic case, one may show that the following Wiener- 
Hopf equations define the optimal conditions (References 2 and 9): 
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\r + WW 

u Q - 76 W = ^ (s) 


r 

r# , v 


1 + — 

_ ww_ 


> 

/y-4- 

ww_ 

ii 

% 

) 


where W = W(s), W = W(-s) and , and have inverse Laplace trans- 

forms which are zero for t > 0. In the event that W(s) is a minimum phase 
system, it suffices to say that (s) , ^ z (s) , and 'f, (s) are analytic in the left- 
half plane. Statements such as ^ 1 tt) - 0 for t > 0 are sufficient conditions 
which, when satisfied, insure a minimum value for the performance index. 

It may occur, for example, that (s) = constant. In this case the sufficient 
condition breaks down and one must resort to a more general condition which 
is both necessary and sufficient (refer to Appendix B). 


In the nondeterministic case, precisely the same sort of analysis 
leads to the following expressions (again refer to Appendix B) 


[r*+E[ww]]u 0 -&£[w} = ^ (s) 


\ 


1 + r£ 


ww 


1 


rie£ ir~=:i * % (*) 


1 fr£ — 

ww 


( ww 
( s ) 


(5-6) 


where ) represents the operation of taking the expected value of the quantity 
in brackets. Again <£‘ [fy ($)] , etc. are required to be zero for t > 0. 

In the multi-control case, if the equations which represent the system 
are given in the first-order form 

% ~ + £u, tj - hly, 

with 00 

Z V ~ J ((j'Qij + u r I2u)dt 

o 

and where F and G have entries which are random variables with the joint 
density p(F, G), the analogous procedure is simply to find the optimal control 
which minimizes the expected value of the performance index as well as the 
expressions for the optimal state variables which minimize the expected 
value of the performance measure. Thus, the sum of the convolutions of 
unknown network impulse responses with the optimal state variables, which 
are set equal to the optimal control, will serve to define the required feed- 
back. The tools for this type of an analysis are already available within the 
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framework of the classical calculus of variations. 


Illustrative Example No. 1 


Experience indicates that the simple first-order system shown in 
Figure 5,4 is adequate to demonstrate the important features of the theory. 



Figure 5.4 First-Order System 


It is assumed that the exact value of a is unknown, being described 
only by the uniform probability density function shown in Figure 5. 5 and 
that %( s) is a step input. 



Figure 5.5 Uniform Probability Density 


In order to specify W , we solve for both the optimal error and the 
optimal control. The Wiener-Hopf equation which must be solved in order to 
determine the error which minimizes the expected value of the performance 
index is 
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Since 


WW 


— = - s z +<x z , 


OO 

(w) -Jj-'* * "')*** 


=■ -s 2 + El{a z } . 



7 


3 



-s 2 + 


7_ 

3 


(- s * f.SZ7)(s+ 1.617) 


The Wiener-Hopf equation becomes, after letting 70 (s) = i/s , 


Z 7 * 

s 


°o-Th 2 ‘T' = 




or 


Let 




(5-7) 


CS4-C. * 

e = 0 i 

0 s(s + cc) 


where C q and are unknown coefficients and substitute into Equation 5-7: 


(- s + *;)((?„ s +(?,)-(- s*f 7 /}) 


= j(s) 


(5-8) 




In order for J o / ^(^)^ 1 W dt - 0 , it is sufficient that the numerator on the left- 
hand side of Equation 5-8 be zero for s = 0 and s = oo . (Refer to Appendix B 
where the case for which ^ ( s) = constant is discussed.) 

At s = 0: oc C 1 =• ^3 

c, -- 7 /^ 

As s — ► <X7 , - C 0 Z must equal -s 2 . 


Therefore, 


4 5 ' 


* In finding the optimal error, a direct method described in Reference 9 is 

being used. 
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The final result is 7 

5 + 

<o OOC 

0 _ 7 7 

S(s + oc) 


Notice that this result gives ot, ) t which is analytic in the left- 

half plane. Hence the given is the optimal one. 


It is important to note that one may employ a root square locus anal- 
ysis (Reference 2) on Equation 5-7 to find the optimal closed-loop roots. This 
would be done for a more complicated example, but it is not at all necessary 
for this simple case since 


oC =■ 




(5-9) 


We proceed now to find the optimal control by solving the Wiener- 
Hopf equation 

(510) 

This will turn out to be a more difficult task in which one is confronted with 
the task of analyzing Laplace transforms which are transcendental functions 
of s. 


The first task is to evaluate the 

/ 


E (ww)= E 


- s z + a 2 


r°° *p(a.)da. C da- 
~~ J *1 -s z + c 


- co 


Thus 




Sf a 


s-c i 


2 1 . 


s 2 -f s- 2 


s*-s-2 


Using the series 


Z 

valid for x > 0 , one obtains 

e { ww ] = - + 


r*-'_ 

lx- 1)3 

i 



t- 

_j 

\x + i) 

5 j 

, v + t) 



1 

+ — 


1 


S z -Z 3 s 2 -2/ 5 l s 2 - 2 


s z g - Z 

S*-s-2 


(5-11) 


(5-12) 


>0 < 5 - 13 > 


Checking the series, it is found that the first two terms afford an excellent 
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approximation to the log function and one may use 


eIwwI - 

1 ' (s>J7) 3 (-s*J7) 3 

Notice that the .fz 1 is the geometric mean of the density function. 

I 


The next task is to evaluate E 

Ep-J. 2 


S + Cl 


l-s + a 


da 


Using Equation 5-12 gives 

/ 


1 


= Mv -s + a ~ JUt 

-5 i- CL 1 


-Sf 2 
- S+/ 


-s + a 


= 2 


- 5 + 


- * (r) 3 f 

I 


' s+ 2 1_ 


M z 


w 


J 


-5*1 


2 *- 


S"35 + 


9* 333 


s*-3 5 + 4 


or 


' ■ 

- ; 

(-s +-jJ 2+ ( ,Z89]Z 

\-s + a 

3 

1 * \ Z 

i J 

S+ 2 

L r'iJ J 


(5-14) 


The Wiener-Hopf equation for the optimal control is now 

(s-yiTf-svrj 3 (-3*1/ 

The root square locus for this system, plotted as | s | vs relative 
damping in Figure 5. 6 shows that for values of 1/ r of approximately one, 
the equation 


(s + 1 . 385) 2 + (. 283) 2 = 0 
defines two of the optimal roots of the system. 

Solving for the optimal control, using the equation (Reference 9) 

a 


u - — 

6 y 


Y 
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where yy a<) - Q . ^ 


one finds, when 1/ r = 68/3, for 


example, 





(-¥)(- s *&) 3 

Jr s *t) + ( 289 ) z 


( 5 + 5.001 ) |j (s + 1. 335)*+ (. Z83)*\ 
4.44 ( s+JT) 3 

S (- s + 5. 00 f)^. s+ -1J 3 (-s + /.38S) z +(.238) i 



s(s + f.OOf)r(s+ t.38S) z + (.283)*' 


Thus 


w c 


since & 0 = 



4-. 44 


s + 5 


S+ 5.001/ 


(W? J 2 | 

(s + 1.385) z + (.Z80) 2 j 


S + 7 /l5 1 

TTI 17 when TT = ~/T~ 
S ( s + 5) ' o 


s + jr _) 

*+■7/5 1 


135 


z s 
DEG 
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225 

1.0 | S , 10 

Figure 5. 6 Root Square Locus for Example 1 
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(5-16) 


(5-17) 


(5-18) 
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If p(a) had the form given in Figure 5. 7, 



the optimal compensating network would have been 

(4.ss)(s + JWs 

(s 4 Vis) < 5 ' 19 ) 


It is apparent that there will be no essential difference between the 
closed-loop system which uses the compensation network given by Equation 
5-18 and the system using the compensation of Equation 5-19. That is, the 
parameter uncertainty was not enough to cause any drastic change in the 
compensating network. 


Illustrative Example No. 2 


It is fairly obvious, from the first example, that the order of the 
compensating network will increase as the information concerning the system 
becomes more '’distributed”. This is not a surprising result and we expect 
it to hold true for all realistic density functions. To illustrate this, the 
optimal compensating network for the density function of Figure 5. 8 and the 
open-loop system of the previous example will be found. 
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The various averages are 

✓ . CK7 


C i ^ 

,^| = / (-s z /-a 2 ) \ $(a-t j + y S(a-z) da = -s 2 + 

E L B L/"J^±i£H* . 

I J - J oo -s z + a. z f-s^f) 


S U A 
5 5 


■s z + a z 


iw] = f [j s{a ^ 4 s ( a ~v\ da 

i J J -s+a 


_ . 6 


1 _oo ~ s + a (rs+t)(-s + 2 ) 

The expression for the optimal control, when ^ = 1/ s, becomes; 


/ e2 8 ) 


i 

( 

r s * tv 

U, . - 

r 

r 


i. J 

(-S z + l)(-5 2 4-4) 


If r = 2/15; — 1 

(s z ^ 4. 265 423s + 4 ) ~y~ (-5*12) _ 

_ (s+l)(s+2)(-s+1)(-5+2) _ 0 sC-St l)(- s+ Z) ^ 

The optimal control is 

Z2s(^^ /X s+ 

15 “ s(s->- 3.333)($ + f.Z) 

The Wiener -Hopf equation for the optimal error is: 

7 .s*(- ***■■%)_ e[-s 2 * 




as + b 

e ° ~ s (s -h 3-3) 


(- 5 +■ 3 . 3 )(as + b)- J- s 2 v- 


*3b- 1 f 

*>' ih 
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2 2 

Everything is well behaved at s = * oo if - as + s - 0. 
Therefore a = 1 


, s + 1.03 

and q - 

° s(s +5.3) 


* 




2. 25 (s + t)(s+ 2 ){'f){s + 3. ?) 

5 + S333)(s + 1.2)(s+ 1.03) 


I 9+1 W 5*.2 W S 5. 3 | 

- Z. 2 S^ tm0S )( 5i . f' 2/(5 + 3.333) 


If p(a) had been as given in Figure 5. 9, 


(5-24) 



Figure 5. 9 Alternate Density for Example 2 


the optimal compensating network would be 


K 



2.22 


(s + 2) 
(s+ 1. 18) 


(5-25) 


which is very similar to that of Equation 5-24. 

The importance of this simple example should not be overlooked. It 
indicates that if one can feel justified in approximating a given density function 
by a weighted set of delta functions, no additional analytical difficulties will 
be experienced in solving for the optimal compensating network than would 
occur in the completely deterministic case. Moreover, it is seen that the 
order of the optimal control expression increases by one for each additional 
delta function which is added to the approximation of the true density function. 
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In addition, all the optimal theory tools such as the root square locus and 
Bode plots (Reference 2), can be brought to bear on the problem, giving the 
engineer a clearer insight as to the effect of any particular uncertain coef- 
ficient on the system’s performance. 


Illustrative Example No. 3 

One intuitively feels that the more realistic density functions for de- 
scribing uncertain parameters in real physical systems should monotonically 
decrease from some peak value and go to zero for two finite values of the 
argument. That is, they should have a form such as that shown in Figure 
5. 10 where the peak value occurs at a = a, and the density function is identi- 
cally zero for a ^ a^ and a — a^. A probable density function which satis- 
fies these requirements and, in addition, does not introduce excessive analy- 
tical difficulties into the problem is the ’’Beta” density function. This func- 
tion is described by Equation 5-26. 



Figure 5.10 A Nonuniform Density 


i> (&) = 


A(a-a 0 ) b (a z -a) c 

0 e^Lse cohere 


*-0 * 


a ± a z 


(5-26) 


The scale factor A is obtained by integrating the density and setting the re- 
sult equal to one. That is, 


Thus 



a Z 

*p(a)da. 


= /. 


P (b + a * Z ) 

r (b+t) ru+o 


(5-27) 


In addition, the peak of the density function occurs for 

ba, + ca 0 

a - — 

' c. + b 


(5-28) 


For illustrative purposes, let b = c = 1, a =1 and a ? = 2. The 

O £ 


71 


IH-2089-F -1 


density function becomes 


-p(&) - &(- a 2 + 3a - 2) 


(5-29) 


Using the system of the previous problem the various expectation 
operations yield (again let fl?(s) = 1/s): 


E l w *h E {^ 7 ^ 


- % l +4 


- s z + 1 


+ (3 s 6j 




; Z -A 5 - 2 


s z - s - 2 


(5 -30) 


The series given in Equation 5-12 yields the following approximation 


for Equation 5-30; 


/ 


-s z + a z 


= W,W f 


where 


(s+ .3S27)\s z + 2(. e )994)( /.SOf)s-h(t.SOl) 2 ] [s 2 f 2(.8W2)(l088)s * (2.0&8) 2 '] 

W = (5-31) 

(s + /T) 3 (s + JW?) 2 


With this information one may follow the analysis procedures used in 
the previous examples to show that 

2.3 

— (5-32) 


s + . 


- 


s/s+oc) 


where 


<C 




results from the solution of the Wiener -Hopf equation 


± 

r 




~ % • 


The equation for the optimal control is now 


1 + - w.w t 




The root square locus for 

f+yWiW, = £ 


(5-33) 
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is given in Figure 5.11. It is seen that when 1/r ^ 10 the equation 
[V* 2(. 9994-)(l.50l)s + (l.50l) Z ~Js Z + Z(.8992)(2.088 )s+(2.468) Z '\ 


U)„ = 1.501 


= 2.088 

% = . 9094 _ 


% - .8992 _ 


defines four of the roots of the optimal control. Rewriting Equation 5-33 

3- S j t — - 

AA ^ ^2 

- LL — = ^ 

' J 


v t v t u * ' 


s D z Tr - 


ior which the solution is 


D , 

U. = 

' A 

one finds 

(s 

II — 

" (JvKp, " 

- sD z A _ 

7 

+ 

7T)> 

s 

00„ = 2.088 
% = . 8992 


u) n = 1.501 
X> * .9994 

(s +.43) (s+s) 


when, for example, 1/ r = 20 


(5-34) 


(5-35) 


Vr N 2 D, ' 


zo(j-)(OT(/W) 3 

L A J 

s -0 

fyilSol) 2 (2.088)* (.43) (S ) 


and 


K * 


io.g(s+JT) 3 (st-feu)* {s + 4.72) 


(s + .43)(s + .48d)(s+5) 


o)„ = 2. 88' 


a >„ = 1. soi ~ 

_ X = .8992_ 


_ % = ,99 94_ 


(5-36) 


since 


when 1/r = 20. 


5+0. $88 
5 (s +■ 4--7Z) 


Thus the Beta density function forces the use of a filter which is 
higher than that found for the uniform probability density function (compare 
Equation 5-36 with Equation 5-18). 
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Figure 5.11 Root Square Locus for Example 3 


Booster With One Bending Mode 


The theoretical development in this section has been restricted to the 
single control variable case. The basic idea behind this approach was to 
keep the theory as simple as possible in the hope of coming to some clear- 
cut conclusion as to whether or not ’’acceptable” control systems could be 
designed using the approach. We have seen, in the previous examples, that 
the theory yields reasonable results which approach the usual ones as the 
density functions become sharper and sharper. Conversely, as the density 
functions become broader, the distributed nature of the information forces 
one to hedge on the design by putting in a higher order compensating network. 
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The single control variable theory is, however, of more than just 
academic interest. One specific problem to which it can be applied directly 
is that of controlling the bending modes in a flexible booster where the only 
control available comes from the gimbaled engines. In this example, trans- 
fer functions obtained from Appendix A will be used. These are listed at an 
appropriate point below. 

The performance index will be 

oo 

W = (5 " 37) 

0 

and will pertain to the block diagram of Figure 5. 12. 



Figure 5.12 Booster Block Diagram 


In Figure 5.12, 0 a / (3 C is the transfer function which relates the 
rigid-body pitch angle ( ) to the gimbal deflection ( (3 C ). 

It is to be assumed that a position gyro is located at x meters and 
that the pitch angle sensed by this gyro is given by the equation 

■ % - % - 1 & 

/ th th 

where. Y- is the slope of the i bending modes and is the i normal 

bending mode. For the purpose of this study, it will be assumed that 


0. ■ 0. - V 


(5-38) 


It is not our intent, in using Equation 5-38, to lightly dismiss the important 
effects of the higher order modes. However, our primary purpose here is 


to examine the nature of the solutions 
selves to only the first bending mode. 

At any rate, 



/3c 


and for this reason we restrict our- 
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The Wiener-Hopf equations which must be solved are; 

q#E{w} = ^ 


(5-39) 


and 





K -*-1, 


(5-40) 


where /3^ = optimal control 
= optimal output 

The compensating network is then 


Kj » 




It will be assumed that everything is known except the value of Y/ , 
which is described only by a probability density function. This is actually 
very close to the truth for a numerical analysis of the equations of motion 
from which the 0 / (3^ and £ ; / transfer functions were derived indicates 
that Y f is the one factor which dominates the uncertainty in the W transfer 
function. At the most forward possible instrument location on the booster 
(x = 1Z2. 5 meters), the data listed in Table 4-1 indicates that 


V= - 


123 ± 25% 


It is now assumed that all the values within this range are equally probable 
and described by the uniform probability density function of Figure 5. 13. 
(Note: the density in Figure 5.13 is actually based on = - . 123 ± 35%). 



Figure 5.13 Probability Density Function for Y/ 


The numbers of interest are 


0C L = -.08 
os | a / 
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and t t 

= — L_ = 11.6 5 

0 ^ 2 . ~ .086 

Other numbers of interest are: 

£ [y/j = -.123 
[e f Yi j] Z = . 01513 
E{(Y/) Z ] =. 0157 

The geometric mean of the density function is J (-■ -38) - .1/524 

Thus, the j^E [ Y t ' J is, to 5% accuracy, equal to E I ( Y f ') 2 J 

As a first task, we elect to solve for the optimal control, (3 C . The 

Wiener -Hopf equation is 

\f + H E [ww}]/3 Co '#<?£ {wj = J 

w = _ Y' — 

A A '' A (5-41) 

LSt A = AL a* JL s A 

(3c Z) D 

For convenience, these transfer functions are listed here (refer to Appendix 

A) 
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we find 
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^ n, - y In, » /v, a^M K 2 ?a/ 2 ^ 

dd 


and cfi-.'p} N ' N ’- £ ( Y f' H N ' N ^ N ' R ‘^ £ { wi’iwz 


DD 


The approximation is made that 

[efX'J]’ ■ e((Y/)‘} 

(N,-E{ Y//A/J {n,-e(Y,'}n,) 


e{ww] = 


D 


p 


or 


£ ( ww J = W t w, 
where A J r E{Y/}N 2 

K = — 


D 


Also, 


£ fw) - £ 


N ( - X' A/ z 


p 


£ fw} = W, 


(5-43) 


(5-44) 

(5-45) 

(5-46) 


(5-47) 


The expression for the optimal control has now become, after division by r 


1+ — W.W, 
r 1 


* "4 - 




(5-48) 


Using the tabulated transfer functions, the expression for is: 


- 2.027 


W, = 




3./ 






2 26M4) 

* 2.3 s + / 


(5-49) 


Prior experience with the bending mode problem indicates that a de- 
sirable situation has been achieved. Note that the complex zeros of the W. 
transfer function have been increased in their natural frequency from the 2 . 3 
rad/sec of the 0//? a transfer function to 3.1 rad/sec. Previous experience 
with this problem indicates that it will now be possible to introduce consider- 
able damping (relatively speaking) into the first bending mode. 


The root square locus for this situ ation is shown in Figure 5. 14*. 

*A gain" we have used a log fs"j vs damping plot similar to those obtained 
from the ESIAC. 
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Equation 5-48 is now rewritten in the form 


for which the solution is 


DP 6 r D 0 




?2(Vr)A/f' 


(5-50) 


If q/ r = 100, from the root locus plot (Figure 5. 14), 

T-(.n) 


IS * 


.013 


"j (tJ 


s\Z 2(.8) 

+ ~ s+ / 


I s ) z 2012) 1 / s ) 

l IT/ * Z.5 SW \3(*0 + 7(5-51) 


The damping factor of the first bending mode has been increased from . 005 
to approximately 0. 12. Ignoring the secondary effects of the root at s = -360, 
the remaining three roots constitute an approximation to a T 1 drift minimum 11 
configuration. (The drift minimum model ideally has a pole at s = 0. ) 

Substituting Equation 5-51 into Equation 5-50, and of course, correctly 
taking into account the scale factor when 12 (s) = 1 / s, one finds 


A.* 


-0*9) 




MsszL,., 

2.3 
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(5-52) 


The above expression is the correct one for the optimal control, even 
though the transfer function has right half plane poles. Refer to Appen- 
dix B for a more thorough discussion of this point. 

Before proceeding to the problem of finding the optimal output {<$^ 1.5 ) 
it is pointed out that the result achieved (within the limitations of the approx- 
imations made) is precisely the same as would be obtained had the value of 
Y^' been exactly known and equal to - .123. This is due to the fact that the 
density function was still sharp enough, and the effect of Y-^ was not critical 
enough, for the uncertainty to enter the analysis in any sensitive manner. 

The next problem is to find the optimal value of 0 12 2.5 which mini- 
mizes the expected value of the performance index. The Wiener-Hopf equa- 
tion is now 
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Since 




n 


= 6.2 x /<?-* 


we will ignore the higher order terms in the series expansion of Equation 
5-57. There is a small element of danger in doing this since N t N z - A/,A/ 2 
has a constant term which is identically zero and therefore the discarded 
terms in the series contain free powers of s. For this reason the approxi- 
mation should be checked as a function of frequency. However, experience 
indicates that the denominator will, as a function of frequency, continue to 
dominate. 


Thus 



DD 

N 1 N 1 - . 123 (N 1 N z -N 1 N z )* .0133 N z N z 
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(5-58) 
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where (^» + is the mean value of Y^' = - . 123. 

If the denominator of Equation 5-58 were expanded, one would see 
that the effect on the coefficients of the term . 0133 N 2 is generally small 

and we may safely set 





(5-59) 


For example, this approximation causes a 1. 5% error in the constant term 
of the denominator of Equation 5-58 . 


Using this approximation, one obtains 


or 



? 

VJp t 



[ 
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At this point, the problem will be terminated, since the Wiener -Hopf equa 
tion has become 


j JL f_ 

$ w,w f 




and the compensating network will be the same as would occur were we deal- 
ing with a completely deterministic system in which Y ^ 123. That is, 

we may as well find the compensating network using the conventional optimal 
theory. 


P 


To summarize the results of this section, 
erformance index 

exO 


the expected value of the 


21 / = 


J[ q lie-0 % )\ r /sf}at 


(5-60) 


has been minimized by finding both the optimal control and the optimal output. 
The slope of the first bending mode was considered to be the only unknown 
coefficient in the equations of motion. The results indicate that one may 
introduce a damping of approximately 12% into the 1st bending mode and the 
required compensation will be essentially the same as that used when a 
completely deterministic system with Y^ 1 equal to the mean value of the 
assumed uniform probability density function. The density, which was as- 
sumed to have a spread of ±35%, was not broad enough to force the use of 
a higher order filter than would normally be used. 


The reader should note that while the example considered a sixth- 
order system, only one uncertain coefficient was involved. However, it is 
not too difficult to visualize that the mathematics would remain essentially 
the same had additional uncertain coefficients been included — namely it 
will continue to use the table of integrals of the rational and irrational 
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algebraic functions to solve the problem. 

Summary 

In this section, simple illustrative examples have been used to obtain 
a "feel" for the types of problems which will be encountered when one applies 
linear optimal control to systems containing uncertain parameters. Towards 
this end, only a single-input and single-output feedback control system was 
considered. 

The approach taken was to assume that the logical thing to do in the 
face of parameter uncertainty was to minimize the expected value of the per- 
formance index. Furthermore, the viewpoint was taken that it made no sense 
to find only the optimal control which accomplishes this minimization since 
one would be forcing an open-loop control law on a system which was original- 
ly postulated to be a closed-loop one. It was reasoned that the best approach 
to the problem was to find the compensating network which minimised the 
expected value of the performance index. However, in the face of the analy- 
tical difficulties encountered in attempting to solve directly for this network, 
the subjective decision was made to solve for the optimal variables at the in- 
put and output of the compensating filter and then define the compensation to 
be the ratio of the two. The indications are that this procedure leads to a 
theory which behaves in a reasonable manner as the information concerning 
the system becomes more distributed. Specifically, the design procedure 
forces one to hedge, by increasing the order of the compensating network, 
as the density functions become broader. 

The basic new features which are added to the linear optimal theory 
in the face of parameter uncertainty are the necessity for the use of the 
tables of integrals for the rational and irrational algebraic functions and the 
problem of approximating irrational functions of s with rational ones in 
order that the root square locus may continue to be of use and that the com- 
pensating networks will be rational functions of s. However, the example 
which considered the weighted delta function density shows clearly that this 
type of density function introduces no additional analytical difficulties into the 
problem. That is, the Wiener-Hopf equations which are involved will, in 
general, be of a higher order than those in the deterministic case but no 
transcendental function of s will be involved. Thus, if the engineer feels 
justified in approximating the given density functions by a weighted set of 
delta functions, the vast body of knowledge concerning the solution of Wiener- 
Hopf equations can be easily modified to solve the uncertain parameter prob- 
lem. 


A specific application of the theory was given in the example which 
considered the design of the compensating network required for a flexible 
booster. It was shown that the uncertain parameter, which was the slope of 
the first bending mode, did not enter into the design in a critical fashion. 
Thus the resultant design would be the same as had one initially fixed the 
value of slope of the first bending mode at the mean value of the probability 
density function. 
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Only uniform, "Beta" and "weighted delta functions" densities were 
considered. It seems important to continue these studies for other types of 
probability density functions. 

It would also be worthwhile to obtain equivalent results for the multi- 
controller problem. As noted in the main text, a good approach to this prob- 
lem would be to find the networks which should be interposed between the 
optimal controls and the optimal state variables. Some preliminary work 
along these lines, not reported here, indicate that reasonable results will 
be obtained. For example, as the density functions become sharper and 
sharper, the ''feedback" networks reduce to the gains that one normally as- 
sociates with the feedback control law u 0 = - . 
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SECTION 6 

THE CONSTRAINED GAIN PROBLEM 

Introduction 

A problem of importance to optimal booster control system design 
is the so-called constrained gain problem. Briefly, for settings of control 
weighting and state- variable excursion weighting, acceptable closed-loop 
system response may be obtained. Once the weighting matrices are selected, 
the optimal feedback gains are uniquely specified. These gains may not be 
changed individually without loss of implementation of the optimal solution. 

Often, some of the feedback gains specified by the optimal solution 
are too large for practical application, in which case, the solution must be 
discarded and a different set of weighting matrices chosen. The new weight- 
ing matrices will cause all feedback gains to change, even though only one or 
perhaps two were too large in the original design, and the new selection may 
still not satisfy the constraints on the gain magnitudes. 

A more desirable approach to the constrained gain problem may be 
to directly constrain those gains that are considered too large (in the optimal 
solution) and then to determine the new optimal solution with these con- 
straints imposed. Although large amplitude constraints on the feedback 
gains might solve the problem, it is unlikely, in that case, that any design 
aid such as the root square locus could be evolved. Moreover, the designer 
might have to resort to a dynamic programming solution with its inherent 
difficulties. Rather than work from the hard-constraint approach, a soft- 
constraint approach was undertaken on this project. In this case, the gains 
are constrained by a quadratic measure, or they are specified a priori with 
the help of the usual optimal regulator solution. The objective is to carry 
the analytical approach as far as possible before resorting to numerical or 
heuristic approaches. 

In the study conducted on this project, a number of analytical approaches 
were worked out in detail. However, even the most promising of these is too 
difficult to use in a practical application, and therefore an experimental or 
heuristic approach was also tried. The major results of both the theoretical 
and the experimental approaches are described in this section. 


Theoretical Development 


Approach 1 

The first approach will allow an analytical solution of the problem, 
which is difficult to use, and the feedback gains are apparently not constant. 
It is described briefly below. 


85 



IH-2089-F-1 


Consider a system of the form 

y, - h'ty 4- Cat , cf - 14 10 , ^ 

where the dimensions of the matrices are 

> (» * i) G (» x f) hi (r k n) y (r x ;J 

F (yi x ri) It(jpxl) ki(^xYl) 

Let 

a * [4; a»i --- 


u, 


a, 1 



1 

^1 

I 

* 

id H 

— T — 

• 


* 

• 


• 



— i 

i 


( 6 - 1 ) 


where each partition quantity has the dimensions 
G q (»xl), hL q (lxl) and ^(lxv>) 

Then 


Gll = G, u, + G z v- z + G_p Uf, 

and 

The performance index to be minimized is 

oo 

2\/ ~ f\^tj'Qy ■> l^,* + ^ 2 S Z id’ + ... + klfSp (6-2) 

0 

where S 1 , S 2 ;••• 5^ are positive definite diagonal ( v * n ) matrices. As 
long as these matrices are all positive definite, the control u is implicitly 
constrained. In addition a constraint on any individual feedback gain in 
may be increased by increasing the corresponding values of the diagonal 
element of the matrix. 

The performance index of Equation 6-2 may be extremized subject 
to the constraints of Equation 6-1 by means of the calculus of variations. The 
Euler-Lagrange equations for the extremals are 

G'Qhix +(F- G 1 ^ - G 2 !d t - . . . - £ + % = 0 

\ + G y % ~ 0 * — /| 2| • • ■ j ^ 


(6-3) 
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where % is an ( » * / ) vector Lagrange multiplier. Rearrangement yields 

A7 'QH% + S^x (Q\ %)\ i --0 

- S^^x } (\ - u z* • •• , 'p 

These equations, together with the equations of 6-1 can be solved for % and 
X and then for each . However, the feedback gains are functions 

of the state variable % so that apparently these gains will be significant- 
ly time-variable. Moreover, the equations of 6-4 are nonlinear, making a 
general solution difficult or impossible. 

This approach to the constrained gain problem has been described 
because it exhibits the difficulty of the problem. A major conclusion that 
can be drawn is that the class of admissible feedback gains must be limited 
to those that are constant (in solution times). If the class is not so limited, 
the feedback gains will vary significantly — a situation that may be undesir- 
able for booster design. It may also be concluded that the feedback gain 
constraint problem leads to the solution of nonlinear equations, whether they 
be algebraic or differential. Therefore, every effort must be made to cast 
the problem in its simplest form so that complexity is minimized. 



Approach 2 

The second approach to the gain constraint problem incorporates in 
the problem statement the specification that the feedback gains are to be 
constant. In this approach, it is assumed that the usual optimal regulator 
problem has been solved, and that some (but not all) of the feedback gains 
are too large in magnitude. The feedback gains that are too large are then 
reduced to values that are at the limit of the magnitude constraints, and are 
thus specified at the outset of the constrained gain optimal regulator solution. 
The problem is to determine the remaining gains (those that did not violate 
the magnitude constraints in the usual optimal solution) so as to minimize 
the original performance measure. 

Again consider a system of the form given in Equation 6-1, but assume 
that the dimensions of the matrices are 

% (m t) Q (n x i) y (r x i) \6 (f * r>) 72 (t * >) 

b (r?*r>) U (] y t) H (f * n) Q (V * r) 

The performance measure to be minimized is 

oo 

2 V - J Q y y- u> f 13 u) dt 

o oq ( 6 - 5 ) 

= /V <- IS'g/J] x dt 

o 

where & is considered to be a vector of constants. Some of the elements 
of /fc' will be specified a priori by virtue of the magnitude constraints. The 
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remainder of the elements must then be determined in such a way as to mini- 
mize the performance measure of Equation 6-5. 

The minimization process is carried out by expressing the performance 
index in terms of Id and then setting equal to zero the expression for the 
partial derivative of the measure with respect to each element of Id that is 
to be determined. The Laplace transform of % it) is 

X [%(iP\ = kL) = Xls-F+G&i 1 ^ ( 6 - 6 ) 

where % Q is the initial condition vector. If use is made of an extended 
form of Parseval’s theorem, the performance measure may be written as 

S'(-s f /6)3(^)T(s, ld)dS (6-7) 



where 


B(/d) s M'QH* WX 


To obtain the extremal values of V for any given gain, the partial der- 
ivative with respect to that gain is taken and set equal to zero; 


dV 

<9£/ 


= 0 


where 'i- is one element of the vector id 

Equation 6-8 may be evaluated by making use of the equation 

— r r 1 s, x) b (*) x ( % A = -L [ Z ' '(- s,ic) B ( A4) X ( s, k)] 




[x7-5,/4)B(^)X(s, *)] 




( 6 - 8 ) 


(6-9) 


where id 0 is equal to Id , but is treated as a constant while the partial 
derivative is being taken. In the derivative operation, it is necessary to 
differentiate inverse matrices with respect to . The operation may be 

performed by using the formula 



where 

Mb M(5, )6) h [l s -r + ££] 


( 6 - 10 ) 


After a great deal of mathematical manipulation, it is found that the condition 
on "i-i for minimization of the performance measure is 
j°° 

-L f el I(S,kZ)X'(-S,U)[p(V)M- 1 (S i k')G-k''ie']c/S = 0 (6-1!) 

J -jo* 
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th 

where 0 Z * is an tf vector with zeros for all elements except for the i 
element, which is 1. 


The optimal condition yields a scalar equation in the frequency domain. 
Equation 6-11 can be rewritten in the form 




ds =0 

( 6 - 12 ) 


where 


\m(s,h)\ is the determinant of M(5 > 1 d) 


and N(S,\t) 


is such that M' 1 (5, It) = 


N(S,jd) 
lM(S, K)\ 


If the closed-loop control system is to be stable, then all of the roots of 
Ms, i d)l = 0 in S must be in the same half plane. Similarly, all of the roots 
of|M(-S,fc)| = 0 must be in the opposite half plane. Moreover, if the left-hand 
side of Equation 6-12 is to be equal to zero, the total integrand must be 
analytic in one or the other of the two half planes. This condition can be met 
if, and only if, the bracketed factor in the integrand of Equation 6-12 contains 
a factor | or two factors |M(S } £)\ . Therefore, the condition 

for the optimal value of the gain *£/ is 

9'. /v(s , kZ) % 0 r' 0 N\~ 5, it) 

under the condition that 


T3(it)N(5, lt)G-k!%\M( SX) 


= 0 


(6-13) 


a. IM(-S,£)I =0 
or 

b. )M(5,lt)l Z = 0 

One equation such as Equation 6-13 can be obtained for each feedback gain 
that is to be optimally adjusted. Since the remaining gains are to be speci- 
fied a priori, the number of conditions obtained is equal to the number of 
optimal gains to be determined. 

Unfortunately, the equations obtained are difficult to solve, except 
for very elementary problems. It is possible that further study of the analy- 
tical solution as given herein may yield additional simplification. Spectral 
factorization, for example, if properly applied might yield a simpler ex- 
pression for the optimal gains. 

Approach 3 

The third theoretical approach makes use of scalar differential equa- 
tions and previously developed tables for optimal feedback gain solution. It 
is, in some ways, a scalar equivalent of Approach 2, but is less general. 
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Consider a single controller system described by the differential 
equation 


u-- ■ a 0 % 


+ 16 +- 


+ a. 


<»> 




(n-2) 

, V (o) - O 


Lm-l) 

% (o) 


given 


(6-14) 


The control variable u is to be a weighted sum of the state variables; that is, 

(6-15) 

Some of the 'i-i 's will be specified a priori because they have exceeded the 
magnitude constraints when the usual optimal regulator solution was obtained. 
The remaining gains are to be chosen so as to be optimal. 


Choose the performance measure: 


* I % * 2 + • •• * <*„./ * l)z + r 0 &,***+•.. + r„. f i** )dt 

O 

OO 

~~ I [(% * r ° *)* 2 + (% + 2 4 • - *(l »-/ * r t,-i *»-, V v ,>2 ] dt 


(6-16) 


where the fa 's and t*i 5 are non-negative constants. 

If use is made of Parseval's theorem, the performance measure may be 
written as . _ 


V = “ J T + >? '£-i Z ) S ** X(V) X^-s) ^5 


(6-17) 


where ^ is the Laplace transform of %(t) ox % . 


The performance measure may be written in the form 


v y - J_ w 

' *„*[£%? w J *.(*) *„(-*) 5 

where ^ J< * 

3 £ (-lyiti+ri *i l )s tl ‘ 

i z 0 

and 

( s ) = (&o+ ^o) + (&1+ £f) 5 + + (&r)-i + ^Y)-f) S ” + s” 


(6-18) 


The performance measure is now in a form that may be evaluated using the 
tables of Appendix D in Reference 8. These tables do not require the factori- 
zation of the numerator and therefore are well suited to this problem. 


Once the performance measure has been evaluated, it may be mini- 
mized by partial differentiation with respect to each free gain Li and each 
result set equal to zero; that is, 

~ibfc7 = ° > 0 “ L ~ (6-19) 


90 


IH-2089-F-1 


Once again, one equation can be obtained for each feedback gain that is to 
be optimally adjusted. Since the remaining gains are to be specified a priori, 
the number of conditions obtained is equal to the number of optimal gains to 
be determined. 

All three of the previous theoretical approaches are very difficult to 
use in practice when the order of the system is greater than three or four. 

The booster design problem, when bending modes are considered, is of 
considerably higher order than third or fourth. Theoretical approaches, such 
as those developed herein cannot be considered as adequate for solving the 
booster constrained gain problem because they are too complex. 


Experimental Solutions and Applications 

To obtain a practical solution of the booster constrained gain problem, 
it is apparent that one must resort to heuristic programming methods. In 
this case, one performs a systematic search for the minimum value of the 
performance measure. Either a digital or an analog computer may be used. 


An experimental study was performed on this project to determine 
the usefulness of a given heuristic method. The dynamics of a booster con- 
sidering actuator and one bending mode were programmed on an analog 
computer. The optimal regulator feedback control system and the computation 
of the performance measure for return from an initial condition were also 
programmed. The equations of motion and performance index used are given 
below: 
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( 6 - 20 ) 

( 6 - 21 ) 


In order that minimization of the performance measure might be accom- 
plished, the computer was switched to the repetitive operation mode. Then 
the feedback gains were adjusted until the final value of the performance 
measure was minimized. Each parameter was adjusted in sequence until 
there was no visible further reduction in the performance measure. 


Initially the computer was programmed for the usual optimal regu- 
lator configuration. A check of this condition on the computer showed the 
analog computer values for the feedback gains to be very close to those com- 
puted by the Ricatti equation solution. After the optimal regulator problem 
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was programmed, one or more of the feedback gains were arbitrarily reduced 
in value and the remaining gains were readjusted so as to minimize the per- 
formance measure. In this way, optimal solutions with constrained gains 
were obtained. The details of the various experiments are described below. 

In the first experiment, the feedback gain from the state variable 
was constrained at various values. In the unconstrained condition, the 
value of feedback gain from is -3. 74. Accordingly, solutions were ob- 

tained for ^ 2 = -3. 74, -3. 00, -2. 50, -2. 00, -1 . 00, and -0. 50. Table 6 . 1 
lists the optimal feedback gains for the various values of constraint on . 

It is seen that there is little effect on the other feedback gains when is 
constrained. However, the locus of closed-loop poles is heavily affected by 
the constraint on % z . This locus is plotted in Figure 6.1 for each value 
of ^ . The locus shows that the major effect of the constraint on the opti- 

mal closed-loop system is to reduce the damping of the resonant pole pair 
associated with the rigid-body dynamics. There is little change in natural 
frequency of this pole pair. Additionally, the damping of the bending mode 
poles is decreased somewhat as becomes heavily constrained. The re- 

sults of this experiment show that analog computer solutions for the optimal 
gains under constraints produce accurate results. There is little deviation 
of the locus of closed-loop poles from the smooth curves of Figure 6 . 1. 

The second experiment involved constraints on two feedback gains. 

First the feedback from the actuator state variable (3 was set to zero, 
thus ^ = 0. Then the feedback from 0% , that is, ^ z , was again con- 

strained at various values. Table 6.2 contains the optimal values of the re- 
maining feedback gains: , and it $ . It can be seen that the 

feedback gains are affected to a slightly greater extent when £4 = 0 than 
when it is optimally chosen. Once again, the locus of closed-loop poles is 
greatly changed as a function of the constraint on )t z , as is seen in Figure 
6.2. Both the damping and the natural frequency of the resonant pole pair 
of the rigid-body dynamics are reduced. It is clear that the lower values 
of ic 2 produce an optimal solution that has low damping of this pole pair. 

In addition, the damping of the bending mode poles is decreased for low values 
of . The locus, again, fits smooth curves accurately. 


TABLE 6.1 


OPTIMAL FEEDBACK GAINS WITH CONSTRAINT IMPOSED ON t z 


RUN NO. 

% 

^ z 

% 

£4 

*5 

% 

1 * 

-1 . 98 

-3. 74 



. 0653 

. 0216 

. 1 31 

2 

- 1.96 

-3. 00 

370 

. 0684 

. 0242 

. 094 

3 

-1 . 84 

-2. 50 

368 

. 0712 

. 0239 

. 083 

4 

-1 . 93 

- 2 . 00 

364 

. 0656 

. 0270 

. 104 

5 

- 2 . 02 

- 1. 00 

369 

. 0663 

. 0191 

. 124 

6 

-1 . 98 

-0. 50 

364 

. 0684 

. 0147 

. 129 


*No constraint on in Run 1. 
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Figure 6,2 Optimal Closed-Loop Pole Locus 
as a Function of the Constraint on = 0) 
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TABLE 6.2 

OPTIMAL FEEDBACK GAINS WITH CONSTRAINT IMPOSED ON , 
AND % SET TO ZERO 


RUN NO. 

% 


^3 


% 

*■, 

7* 

TZ753 

— rrn — 

376 

. 0658 

. 0203 

0 

8 

-2. 33 

-3. 00 

369 

. 0696 

. 0208 

0 

9 

-2. 20 

-2. 50 

369 

. 0708 

. 0230 

0 

10 

-2. 10 

-2. 00 

367 

. 0754 

. 0249 

0 

11 

-1.99 

-1 . 00 

-. 361 

. 0755 

. 0224 

0 


*No constraint on in Run 7. 

The third experiment was similar to the first, except that was 
constrained instead of ^ 2 , . It was noted in obtaining the previous solutions 

that the value of the performance measure was very sensitive to the setting 
of i: 3 , which is the feedback from the state variable, CC . Therefore, one 
would expect large changes in the gains and in the closed-loop locus when 
is constrained. These large changes do, indeed, occur as Table 6.3 
and Figure 6.3 show. All five remaining gains undergo large changes as 

is constrained. The value of ^ changes sign. The resonant pole pair 
associated with the rigid body has a greatly reduced natural frequency, al- 
though the damping is not changed significantly. Moreover for the bending 
mode poles, the damping is reduced slightly and the natural frequency is in- 
creased slightly. It would have been desirable to obtain more solutions for 
larger constraints on ^ 3 , but this would have required complete rescaling 

of the problem on the analog computer. 


TABLE 6. 3 

OPTIMAL FEEDBACK GAINS WITH CONSTRAINT IMPOSED ON 


RUN NO. 

% 

^ z 

*3 

u 

*5 


1 * 

-1 . 98 

-3. 74 

-. 370 

. 0653 

. 0216 

. 131 

12 

-0. 49 

- 1 . 66 

300 

. 0643 

. 0158 

. 112 

1 3 

-0. 37 

-1. 19 

200 

. 0160 

. 0115 

-.017 


*No constraint on ^ in Run 1. 


Conclusions 


The studies of the constrained gain problem that have been performed 
on this project have made certain aspects of the problem clear. First it can 
be seen from the complexities of the mathematical approaches that analytical 
solutions are of little value at this time. If the problem is expressed in its 
simplest form, a set of nonlinear, simultaneous algebraic equations must be 
solved. Since there are no general methods for solution of these nonlinear 
algebraic equations, the analytical solutions are too difficult to allow straight- 
forward solution. It is possible that certain specific solutions could be worked 
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out by substitution and tedious computation, but these methods are not recom- 
mended for high order systems. 

In order that heuristic programming methods might be evaluated in 
regard to application to the constrained gain problem, a sample booster prob- 
lem was programmed on the analog computer. The performance index was 
minimized without difficulty when various gains were constrained. When the 
usual optimal regulator solution was used as a starting point, and one or 
more of the gains were gradually reduced so as to approach a desired con- 
straint, a minimum could be found. No problems were encountered that 
could not be traced to the computer itself or its scalings. 

It is clear that in some cases no constrained gain solution exists. 

For example, if tlm upen-loop system is unstable, it may be possible to con- 
strain one or more of the feedback gains (in magnitude) in such a way that 
there is no stable closed-loop solution. Physical reasoning will generally 
make the designer aware of such situations, and of course, they must be 
avoided in any heuristic programming method. 

The experimental study has demonstrated the feasibility of developing 
a digital computer program which automatically performs a systematic 
search for optimal constrained gain solutions. A strategy could be developed 
for adjusting each parameter until the performance measure is minimized. 

If the program uses initial gain settings taken from the usual optimal regula- 
tor problem and then gradually shifts the constrained gains, an optimal solu- 
tion can be obtained for most practical cases. The program could be developed 
in such a way as to accept matrix equations of motion and arbitrary quadratic 
performance measures similar to those used in the usual regulator problem. 

It is recommended that such a program be developed, since it may solve 
most constrained gain problems effectively and easily. 


I 
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SECTION 7 

CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 

1. Linear optimal control and associated matrix techniques can be used 
in a systematic and effective way to conceptually design satisfactory control 
systems for large, flexible launch vehicles. 

2. The use of models in the design procedure simplifies the selection of 
suitable quadratic performance indices. 

3. The optimal feedback control law for the flexible launch vehicle can 
be obtained in a direct, straightforward manner, not requiring the solution 
of a matrix Riccati or Wiener -Hopf equation. 

4. The control law can be expressed in terms of quantities sensed 

on the body of the vehicle. The selection of the sensors and their locations 
is an important design consideration, for the feedback control law may not 
be realizable in all cases. 

5. The optimally controlled system appears to be insensitive to variations 
of the bending mode shapes and slopes. 

6. Optimal compensation networks can be specified for systems with 
parameters whose values are known only as random variables described on 
a sample space. 

7. The ’’constrained gain” problem is difficult to solve in an analytical 
manner. However, it appears that solutions can be experimentally obtained 
quickly and easily using an analog or digital computer. 


R ecommendations 


1. It has been demonstrated that linear optimal control techniques can 

be effectively used as a design tool to conceptually specify satisfactory feed- 
back control systems for flexible launch vehicles. The techniques are rela- 
tively new, however, and all the unique aspects of the methods have not yet 
been investigated. In essence, the wealth of experience and background in- 
formation collected for conventional design methods over the past three 
decades is not directly available to the optimal control system designer. A 
more measured and deliberate pace is therefore required when using optimal 
design methods. However, because of the promise and progress that these 
techniques have shown, optimal control studies for application to the flexible 
launch vehicle control problem should be continued. 

2. Optimal control law synthesis problems should be investigated in more 
detail. Other control law approximation techniques should be studied with 
the objective of defining the extent to which approximations can be made 
without causing an unstable configuration. 
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3. The results concerning parameter variations are promising. A 
more thorough investigation should be conducted to extend the theory of opti- 
mal design of systems subject to parameter variations. 

4. One of the important results of this study is the experimental deter- 
mination that the optimal system is relatively insensitive to variations of 
the bending mode shapes and slopes. These results should be placed on a 
more firm theoretical basis for the elastic launch vehicle. 

5. An investigation of the constrained gain problem should be continued. 
Relationships among the parameters in the performance index and the magni- 
tudes of the feedback gains should be established. Relationships that exist 
among the closed-loop dynamics of an optimal system with one or more feed- 
back gains fixed and with a preselected performance index should be established. 
A digital computer program may be the solution to this particular problem. 
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APPENDIX A 

EQUATIONS OF MOTION AND TRANSFER FUNCTIONS 


The equations of motion used throughout this report were extracted 
from a Marshall Space Flight Center memo labeled "Model Vehicle No. 2 
For Advanced Control Studies". The portions of the memo that are used in 
the body of this report are included in this appendix. The simplifications 
that were made are described and the pertinent transfer functions are in- 
cluded in this appendix. 

The rigid-body coordinate system and the first elastic model geo- 
metry are shown in Figures A. 1 and A. 2, which are reproduced from Refer- 
ence 1. Using the nomenclature of these two figures, the equations of motion 
are 
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Flight Path Curvature Equation 
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Figure A. 2 First Bending Mode Geometry 
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The engine-actuator dynamics satisfy the differential equation 

f3 + 23. Ij6 * 2 259 f3 * 3t, /JO f3 - 31, 130 /3 C 

In addition to the equations of motion given above, the following 
subsidiary equations, which define the quantities measured by an attitude 
gyro, a rate gyro and accelerometer are used. 


Output of an accelerometer 
Rate Gyro 
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Attitude Gyro 



The entire analysis in the report is performed at a fixed operating 
point during the boosters’ launch trajectory. The point chosen was t = 80 
seconds along the nominal trajectory. This operating point is within two 
seconds of the maximum dynamic pressure condition. 

A number of simplifications and approximations were made to elim- 
inate the minor terms in the equations of motion, when these minor terms 
add little in the frequency range of interest. The simplifications are; 

1. The actuator and engine dynamics were approximated by the first- 
order equation 

(3 + 17 . 9(3 - 17.9 ( 3 ^ 

The constant ijT - 17. 9 represents a first-order approximation to the fre- 
quency response and phase characteristics of the third-order engine dynam- 
ical equation. 

2. The inertia of the gimbaled engines is small compared to the total 
inertia of the vehicle, therefore ^ 

3. The pitching acceleration due to engine reaction force is small com- 
pared to pitch acceleration due to control forces, i. e. , 

V'-lcg r- x s £ 4 F- X s £ 

— >> 4 therefore r — 0. 

m Iy ’ rn Iy 

4. The incremental pitching acceleration and normal acceleration 
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caused by vehicle flexibility are negligible. Therefor 
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Using the data given in the Model Vehicle No. 2 memo and using the 
assumptions given above, the equations of motion of the launch vehicle at 
t = 80 seconds were found to be 


(fyr .07330. + .45 (3=0 

.0405 <fi z + oo + . 0)067 X + . 021060 = 0 
(3 + 17,9 0 = 17.9 01 

-5.4S3 0C + + .02317%, + 5.37%, - 15.83 0 = O 

-2.36a, * % 2 * .0564 % z i 31.8 % 2 - 22. 770 = 0 
- 11.8ol + r[ 3 + .0918 j 3 + 84.25 % 3 - 26.250 = 0 


pitching acceleration equation 
flight path curvature 
actuator dynamics 
1st bending mode 
2nd bending mode 
3rd bending mode 


Most of the conceptual design work reported upon in the body of this 
report included two bending modes in the equations of motion. During the 
course of the design work, certain transfer functions were required. These 
transfer functions, which include only two bending modes, are given below; 
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APPENDIX B 

THE WIENER -HOPF EQUATIONS 


This appendix discusses the derivations which are pertinent to the 
theory contained in Section 5 of the report. First, a brief review of the 
derivation of the Wiener-Hopf equation, which specifies the optimal control 
will be given for the case where the plant is precisely known. In the deriva- 
tion, the effect of having a plant with non-minimum phase components will 
be emphasized and then an example will be given to demonstrate the various 
aspects of the problem. After this, the Wiener-Hopf equation which defines 
the optimal control that minimizes the expected value of the performance 
index will be derived. 

It seems necessary to carry out the proofs for only the optimal con- 
trol. Those readers who are interested in deriving, for example, the ex- 
pression which defines the optimal error may mimic the procedure used 
for finding the control. The one difficult problem which seems to affect the 
optimal error theory to a greater extent than it does, for example, the opti- 
mal control theory, is the point of whether or not the sufficient condition for 
a minimum is satisfied when 'fy(s) is a constant. This point is considered at 
an appropriate point in the development. 

The Optimal Control in the Deterministic Case 


Consider the block diagram of Figure B.l 
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Figure B. 1 Closed-Loop System 
In Figure B.l, W represents the fixed elements of the system. 
The performance index 

2V = / (e z +ru 2 )dt , 


(B-l) 


in which r is the positive weighting on the control, is to be minimized by 
solving for the optimal control. 


form 


since 


Equation B-l can be rewritten in the equivalent frequency domain 




2 \/ = 


2/r/ 


*?- Wu)(ie- Wu)+ruu\ds 
y(s) = W(s)u(s) 


(B-2) 


(B-3 ) 
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In Equation B-2, W = W(s), and W = W(-s), etc. 

One now takes a variation on Equation B-2 by letting 

Lb = U 0 + %U, (B-4) 

where Li a is the optimal control and U 1 is physically realizable, but other- 
wise arbitrary. (If this derivation were carried out in the time domain, it 
would probably be more evident that the minimuzation procedure also requires 
that (1, Umc) = a 1 (tmeo) = 0.) The result is 

2V= J a + %(J b + + (B-5) 

where 

J oO 

<4 = ~zWj I I ft- -W* c ) + ra 0 zA ds 

J -j 00 J 

7 b = -^f \(r + V/W)u 0 ~WK}u,ds 

j 0 * 

Jc= ! Jr f {(r*W9)«>-V*U*s 

•J -joa ' J 

j°° 

J 


J a is the optimum value of the performance index, > 0 and = J a since 
J b (-s) = J & (s). From this it follows that a necessary and sufficient condition 
for ti 0 to give the lowest value of J is, for arbitrary u f , J c - 0 (see, for 
example. Reference 4). 

Thus one must consider the equation 

joo 

(r f WW)u, d - WfsJ u f ds (B-6) 

in which all the poles of U i are in the right-half plane. It is at this point that 
care must be exercised because one is tempted to say that a sufficient condition 
for Equation B-6 being identically zero is that all of the poles of the expression 

(r+ \Nvi)u 0 -vne 

must also be inside the right-half plane when the path of integration is com- 
pleted to the left. This sufficient condition is usually expressed as 
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(r*WW)u g -W1S»^(s) (b-7) 

where /j.(s) is analytic in the left-half plane and, moreover, it assumes that 
the integrand of Equation B-6 has a sufficient excess of poles over zeros such 
that the integrand is forced to zero along the left-half plane contour. Equa- 
tion B-7 is a scalar Wiener-Hopf equation for which one will obtain the correct 
answer by requiring that 'f-(s) have right-half plane poles only when W(s) is 
minimum phase. In the event that W(s) is nonminimum phase the correct 
sufficient condition is 


X~Xy(s)]- 'yW-o Ut±o 


for then 



(B-8) 


(B-9) 


is identically zero. The basic difficulty is due to the fact that, when there are 
poles in the right-half plane, W(s) has an impulse response which goes to oo 
as t oo while the impulse response associated with W(-s) goes to oo as 
t — - oo . The danger exists that one will attempt to solve the problem in a 
manner that identifies all right-half plane poles with time functions which 
exist only for t < 0. Clearly a right-half plane pole can also represent an un- 
stable system for t > 0 and one must be careful to associate poles with the 
correct type of time function. Taking these precautions, one may now rewrite 
Equation B-7 in the form 


\r DD + Nn] 
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(B - 1 0) 


where N 

W “ — . 

D 

D represents the open-loop poles of the system and will contain right-half 
plane poles when the open-loop system is unstable. Nevertheless, the time 
function associated with D is identically zero for t < 0 while the time function 
associated with D is zero for t > 0. [rDV +■ NN ] is now written as the product 
of a right-half plane component and a left-half plane component, that is, 

[ rVD+NN ] = (r DD +-nn) + (rVD +Nbi)~ = A A 

where, for example, A gives rise to a stable time function which exists for 
t > 0. That is, we insist that the poles be in the left-hand plane. 


Now define 
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[Nig/ a] has been decomposed into a sum in which the first term on the right 
is associated with a "positive time" function and the second term is associated 
with a "negative time" function. Moreover, because of the definition of A > 
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will be a stable positive time function. In making this statement, 
one assumes that ^ has left-half plane poles. This decomposition can be 
achieved using a partial fraction expansion. 


One can now set 


“o = 
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Nt2 


(rVD +/V/v) + L ( r D D f a/a/9'J 


(B - 1 1 ) 


and verify that this is truly the optimal control by substituting back into Equa- 
tion B-10 to obtain 
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This is a function for which 


Equation B-ll is correct even when W(s) has poles and zeros in the 
right-half plane and is also correct when W(s) has a time delay if one inter- 
prets it properly. That is, 
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(B-13) 


(We will avoid the trap of talking about the "prediction" problem which is 
sometimes incorrectly stated as 


W(s) = — & sT 


since 


X' f \W( sj] ^0 fir t ±0). 


To illustrate the preceding discussion, consider a simple example for which 

(B-14) 


(-s+1)e sT 1 

W( s) = — and i2 - — 

s(-s-bZ) 5 


Let N = (-s+)) and D= S (-?■<■ z)e sT and use Equation B-l 1 ; 

s(-s + 2)e sr 
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where 


and 


b = 


SF 



1 


The only "left-half plane pole" of concern in the partial fraction expansion is 
due to the step input. Therefore 


where 
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If one were to use the 
blindly, the result 


"analytic in the left-half plane" 

S-h Z 

Jr (s z fd5 + b) 


requirement 


would be obtained. This does not satisfy the Wiener-Hopf condition. 


Substituting Equation B-15 into the Wiener-Hopf equation 


gives 
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(B-16) 


Thus f. , ['^(s)] = 0 for t ± 0 b ecause ( S + 1 ) G sT is, by definition, a iactor of 
D and gives rise to a time function which exists only for negative time. 


In any situation where there is some doubt as to what to do (for ex- 
ample, doe s one as sociate e ~ sT with N or e 4 sT with D ? ) the only recourse is 
to check the possible alternatives and see whether or not Equation B-9 has 
been satisfied. 

It may occur that ^ (s) turns out to be a constant (K) which is, of 
course, analytic in the entire s plane. In this event, the sufficient condition, 
at least in its more general interpretation, breaks down, because 
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00 


— - — J /}(s) u 1 (s) 0/5 = J ^ (t) u 1 (t)dt 
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= lcf Js(t)u,(t)dt ~u t(ts 


o) 


It is thus necessary to evoke the additional requirement that tl. -u.. , = 0. 

The results for the single variable case are easily extended to the 
multi-control case {see, for example, Reference 2) when nonminimum phase 
components are involved. 

The Optimal Control in the Random Variable Case 


In the random variable case we again start with the performance index 

OO 

ZV= J ru z )dt (B-17) 


which is also equal to, by Parseval’s theorem, 
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Therefore ^ 

£ f 2\/l = • E 

1 J 

Since E { 2V ] really means that one must multiply 2 V by the multi-dimensional 
probability density function associated with all of the unknown system param- 
eters, and integrate over the proper domain, the conditions under which one 
can interchange the order of integration and the expectation operation reduce 
down to the sufficient conditions one considers for inverting orders of inte- 
gration in iterated integrals. For example, the formula 




is valid if both sides exist and at least one converges absolutely. That is, at 
least one of the integrals 
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is finite. 


Thus if we wish to write 
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where P(A) represents the multi-dimensional cumulative probability function 
for the unknown parameters (described by the vector A), must exist (i.e., 
be less than infinity). 

The writers know of only one practical case for which this condition 
is not satisfied. If the input to a system is such that a finite following error 
must result, then Equation B-17 is infinite (e.g., a ramp input into a type one 
system). However, even in this case it can be shown that the integrand of the 
performance index has a minimum steady state value. This point is discussed 
at some length in Section 6 of Reference 2. 


We now assume that the sufficient conditions have been satisfied and 

write 
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+ ru,u | ds 
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Assuming that IS 
the expected value of the 


is deterministic and that the control which minimizes 
performance index is to be found, one obtains 
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Again, let 
optimal condition 


U - U- 0 + t substitute and clear through to obtain the 

\r + E[wwj ^ 0 - K?e[wJ ■= (s) (B - 1 9) 


where 


#"'[><*)] “ ftt) ~ 0 -for t ±0 


It is easy to see that interchanging the expectation and integration op- 
erationwill give the same answers as were obtained in the deterministic case 
with WW replaced by E j WW ), etc. Thus the optimal error condition becomes 
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(B-20) 


and the optimal output condition becomes 


(B — 2 1 ) 


Notice that the basic philosophy emphasized in this very straight- 
forward and simple derivation is that W is fixed for any given system but 
not precisely known. We then subjectively decided to take the usual per- 
formance index and find its expected value. It is only fair to point out that 
this idea would lead to trivial results for the case where the system was 
deterministic but the input was a random variable for then 


£ 



(B-22) 


would give the result that U 0 = 0 when IS has a zero mean. In this situa- 
tion, a different philosophy is adopted which uses the performance index 


J f = e z ± ra z 


and asks that one find the minimum of a time average. That is, 

-r 

When this problem is explored, one finds that the reasonable physi- 
cal entity to minimize is the system impulse response [ (£J ] . The vari- 

ous manipulations are then carried through completely in the time domain 
and yield the time Wiener-Hopf equation (see, for example, Reference 1). 
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/ - 0 M (Z,) = 0 (B-23 

0 

hr r±o 

where 

0 C e = autocorrelation function of the input 

= autocorrelation function of the input to the fixed 
elements of the system 

F (T t ) ~ the system impulse response which minimizes 
the E { J } . 

Transforming into the frequency domain one obtains the frequency do 
main Wiener -Hopf equation 

tiuWF . l *)-# < B - 24 

where 

ot 1 [fyis)] = 0 for t > 0 

00 = the power spectral density of the input signal 
= the optimal transfer function. 




